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Preface 


The  purpose  of  this  stud/  was  to  investigate  the  geometric  performance  of  a  modification  to 
the  proposed  baseline  constellation  for  the  NAVSTAR  Global  Positioning  System  ( GPS).  To  assist  in 
the  evaluation,  a  computer  program  developed  by  the  Rand  Corporation  ( 1 7)  was  obtained  and 
modified  for  use  in  analyzing  the  geometric  performance  of  this  pseudoranging  navigation  satellite 
system.  The  geometric  perfxmance  of  the  currently  proposed  baseline  constellation  was  analyzed 
using  this  tool  and  was  compared  with  the  performance  of  the  modified  constellation.  Although 
additional  factors  still  need  to  be  considered  in  selecting  the  optimum  satellite  constellation,  the 
modification  presented  here  greatly  reduces  the  percentage  of  system  outages  that  occur  due  to  the 
poor  geometry  of  available  satellites  and,  at  least  on  the  surface,  appears  to  provide  a  significant 
improvement  over  the  present  proposal.  Perhaps  more  importantly,  it  suggests  that  the 
expectation  of  Walker  (22.-4)  that  'continuous  whole-Earth  coverage  would  be  provided  most 
effectively  by  a  system  in  which  the  distribution  of  satellites  over  the  Earth's  surface  was 
maintained  as  uniform  as  possible,"  is  not  necessarily  true .  Since  practically  every  analysis 
to  date  has  assumed  Walker  s  "expectation"  as  a  basis  for  comparing  constellations,  the  results  of 
this  study  could  prove  to  be  of  significant  value  in  future  analyses. 
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and  to  Professor  C.R.  Edstrom  (AF1T  Math  Department),  for  his  assistance  in  understanding  the 
mathematical  concepts  involved. 
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Abstract 

This  investigation  determined  the  effect  of  changes  in  eccentricity  to  the  orbits  of  the 
proposed  Global  Positioning  System  (GPS)  18-satellite  baseline  constellation  by  analyzing  the 
geometric  performance  obtained.  The  effect  of  satellite  losses  upon  global  coverage  was  also 
examined  with  an  emphasis  on  determining  which  combination  of  remaining  satellites  provided 
the  best  and  worst  cases.  The  potential  of  GPS  for  navigation  of  the  space-based  user  was  explored 
by  analyzing  the  geometric  performance  obtained  for  a  variety  of  user  trajectories  and  GPS 
antenna  beemwidths.  A  computer  program  which  analyzes  the  many  aspects  of  the  geometric 
performance  of  pseudoranging  navigation  satellite  systems  was  used  for  the  analysis. 

The  results  of  this  analysis  Indicate  that  a  simple  modification  to  the  baseline  constellation 
could  reduce  system  outages  on  8  global  basis  by  nearly  508.  The  modification  consists  of 
changing  the  shape  of  the  GPS  circular  satellite  orbits  to  slightly  elliptical  ones,  resulting  in  more 
favorable  satellite  geometry  and  fewer  outages  to  the  user  on  8  global  average.  Further 
consideration  to  determine  Its  feasibility  was  recommended.  The  degradation  of  coverage  due  to 
satellite  losses  was  found  to  be  largely  dependent  on  the  combination  of  the  remaining  satellites, 
and  suggests  that  rephasing  the  remaining  satellites  could  significantly  improve  the  degraded 
performance.  The  potential  for  conventional  use  of  GPS  for  navigation  in  space  was  shown  to  exist 
for  the  low  altitude  user,  but  will  be  very  limited  for  the  higher  altitude  user  due  to  the  present 
GPS  antenna  design.  Increasing  the  designed  antenna  beamwidth  was  shown  to  significantly 
improve  performance  for  the  high  altitude  user.  It  wss  recommended  that  this  modification  be 
considered  in  future  GPS  antenna  design,  if  conventional  GPS  navigation  is  to  be  desired  for  the 


high  altitude  space  user. 


GLOBAL  POSITIONING  SYSTEM  -  A  MODIFICATION 


TO  THE  BASELINE  SATELLITE  CONSTELLATION 
FOR  IMPROVED  GEOMETRIC  PERFORMANCE 

I.  Introduction 

The  Navster  Global  Positioning  System  (OPS)  is  a  space-based,  pseudor8nging  navigation 
satellite  system  that  will  provide  worldwide,  nearly  continuous,  three  dimensional  position, 
velocity,  and  coordinated  universal  time  to  the  suitably  equipped  user  ( 2:226).  Although  designed 
primarily  for  global  navigation  of  a  terrestrial  or  near-earth  user,  the  potential  exists  for 
expanding  its  use  to  the  space- based  user.  This  study  analyzes  the  geometric  performance  of  the 
proposed  baseline  orbital  constellation  not  only  for  the  earth-based  user ,  but  for  the  space-based 
user  as  well.  In  addition,  the  effect  on  geometric  performance  of  modifying  antenna  beamwidth  is 
examined  for  a  variety  of  space-based  user  trajectories,  including  both  low  and  high  altitude 
orbits  ,  intermediate  transfer  orbits,  and  typical  ICBM  trajectories.  The  effect  of  satellite  losses 
on  geometric  performance  and  system  accuracy  is  examined  with  a  special  emphasis  on  the  'best’ 
and  ’  worst*  case  scenarios. 

The  position  accuracy  available  from  GPS  can  be  divided  into  two  multiplicative  factors, 
position  dilution  of  precision  ( POOP )  and  other  "system"  errors  ( 1 7: 1 ).  The  "system"  errors 
depend  on  the  accuracy  of  the  ephemeris  data  and  transmitted  time  from  the  satellites,  ionispheric 
and  atmospheric  effects,  and  various  other  error  sources  as  indicated  in  Table  I.  Since  the  POOP 
factors  depend  predominantly  on  the  user/navigation  satellite  geometries,  they  can  be  analyzed 


independently  of  system  errors,  which  depend  on  a  multitude  of  factors.  This  characteristic  allows 
separate  analyses  of  alternative  orbital  configurations,  user  motion,  and  satellite  losses  sustained 
for  the  purposes  of  comparison  and  choosing  the  optimal  constellation  ( 1 7:3). 


TABLE  I 

GPS  Error  Budget  (2:228) 


ERROR  CONTRIBUTION 
(METERS) 


Space  &  Control 
Segments 


User 

Equipment 


ERROR  SOURCE 


L,/L2 

P-Code 


Clock  &  Navigation  Subsystem  Stability  2.1(1  -Sigma) 

Predictability  of  Satellite  Perturbation  1 .0 

Ephemeris  Prediction  8>  Model  Implementation  2.5 
Other  0.7 


Ionospheric  Delay  Compensation 
Tropospheric  Delay  Compensation 
Receiver  Noise  &  Resolution 
Multipath 
Other 


Ranging  Error  -  RSS  Total  1-Sigma 


3D  RMS  Navigation  Error  -  (POOP  RMS)  (Ranging  Error) 
*  (3.7)  (Ranging  Error) 


Spherical  Error  Probable  *  (0.8)  (30  W1S) 


After  fully  analyzing  the  geometric  performance  of  the  proposed  baseline  constellation,  a 
modification  to  this  constellation  will  be  incorporated  and  then  be  analyzed  for  comparison.  This 
modification  will  "target"  the  geographic  areas  of  weaker  coverage  in  an  attempt  to  reduce  the 


number  and  length  of  system  outages  that  occur  (due  to  poor  geometry)  on  a  global  basis  by 


2 


changing  the  eccentricity  of  the  six  orbital  planes  and  selectively  positioning  the  periapsis  point  of 
each  orbit.  As  the  primary  tool  for  the  analysis,  a  computer  program  on  the  geometric 
performance  of  pseudoranging  navigation  satellite  systems,  developed  by  the  Rand  Corporation 
( 17),  was  obtained  and  modified  for  this  purpose.  (For  more  information  on  the  computer 
program  and  its  operation,  see  Appendix  B). 

In  order  to  better  understand  the  significance  of  the  geometric  performance  of  each 
constellation  and  its  effect  on  system  accuracy,  a  thorough  review  of  the  terminology  and  OPS 
concept  of  operation  is  necessary.  Therefore,  Chapter  II  will  be  devoted  to  providing  the 
background  necessary  for  understanding  the  system  concepts  involved.  The  mathematical 
derivations  for  the  navigation  equations  will  be  presented  in  Chapter  III,  and  the  dilutions  of 
precision  ( DOPs)  will  be  defined  and  related  to  overall  system  accuracy. 

Chapters  IY  and  V  will  be  devoted  to  the  data  analysis  for  both  the  proposed  baseline 
constellation  and  the  modification  to  that  constellation.  The  geometric  performance  of  each 
constellation  will  be  directly  compared  with  that  of  the  other,  and  from  this  comparison, 
conclusions  and  recommendations  will  be  made  and  presented  in  Chapter  VI . 

To  further  assist  the  reader  in  understanding  the  information  contained  in  this  report,  a 
glossary  of  technical  terms  is  provided  in  Appendix  A.  Appendix  B  provides  an  explanation  of  the 
variables  and  subroutines  of  the  computer  program  utilized  in  the  analysis  and  explains  the 
modifications  that  were  incorporated.  The  complete  computer  program  listing  is  provided  in 
Appendix  C,  followed  by  samples  of  the  program  output  (Appendix  D)  and  selected  data  extracts 
(Appendix  E)  from  the  more  than  1 50  computer  runs  made  during  this  evaluation. 


ii.  System  Concept 


Background 

Since  the  early  1 960's,  the  idee  that  navigation  and  positioning  could  be  accomplished  using 
radio  signals  transmitted  from  satellites  has  been  actively  pursued  by  both  the  Navy  and  the  Air 
Force.  Each  service  separately  established  its  own  concept  of  such  a  system  through  an  extensive 
program  of  studies  and  tests  designed  to  demonstrate  the  feasibility  of  a  space-based  positioning 
and  navigation  system  ( 1 9:2 1 . 1 . 1 ).  The  success  of  the  Navy  Navigation  Satellite  System ,  better 
known  as  TRANSIT ,  stimulated  both  the  Navy  and  Air  Force  to  develop  a  more  advanced  system  that 
would  provide  enhanced  capabilities  and  global  coverage.  The  Navy's  concept,  T I  MAT  ION,  was 
essentially  a  two-dimensional  system  aid  could  not  provide  position  updates  in  a  high-dynamic 
aircraft  environment.  The  Air  Force  concept,  Program  621B,  could  provide  the  high-dynamic 
capability  but  had  its  own  shortcomings  os  well,  particularly  from  a  survivability  standpoint 
(10:1177). 

Recognizing  the  need  to  integrate  these  systems,  the  Deputy  Secretary  of  Defense  issued  a 
memorandum  in  1 973  which  designated  the  Air  Force  as  the  primary  agency  to  develop,  test,  and 
deploy  a  single  system  that  could  best  serve  the  needs  of  the  defense  ( 1 9:2 1.1.1).  A  system 
concept  which  combined  the  best  features  of  both  programs,  designated  the  Navstar  Global 
Positioning  System ,  was  the  resulting  system  design.  Management  of  the  Navstar  OPS  was  assigned 
to  the  Joint  Program  Office  (JPO)  located  at  Space  Division  in  Los  Angeles. 

When  first  conceived,  the  design  for  the  fully  operational  system  consisted  of  a  total  of  24 
satellites,  deployed  with  eight  satellites  uniformly  distributed  in  each  of  three  orbital  planes, 
providing  continuous  three-dimensional  global  coverage  with  predicted  accuracies  in  the  ten 
meter  range  ( 1 0: 1 1 77).  In  1 978  aid  1 979  Defense  budgetary  constraints  forced  a  reduction  in 


funding  for  the  GPS  program.  As  a  result,  in  1 980,  the  Navstar  GPS  program  was  restructured 
and  the  number  of  satellites  for  the  fully  operational  system  was  reduced  from  the  24  originally 
planned  to  1 6.  Exhaustive  studies  have  been  made  since  that  time  to  establish  the  optimum  orbital 
configuration  for  these  remaining  18  satellites.  System  accuracy,  survivability,  satellite 
visibility,  ease  of  buildup,  location  and  duration  of  outages,  ease  of  sparing  and  replacement,  and 
growth  potential  of  the  constellation  to  24  satellites  were  considered  in  the  selection  of  what  is 
now  the  baseline  constellation,  which  will  consist  of  18  satellites,  uniformly  distributed  in  six 
orbital  planes  with  three  satellites  per  plane,  and  will  provide  nearly  continuous  world-wide 
coverage  that  optimizes  accuracy  over  the  primary  areas  of  interest  ( 1 5:  E9.3. 1 ). 


Figure  2.  Nevstar  GPS  Segments  (6:2 18) 


Space  Segment.  Designed  to  be  fully  operational  in  late  1 988 ,  the  xbital  constellation 
will  consist  of  1 8  satellites  deployed  in  circular,  1 2  hour  orbits  of  26,600  km  radius  ( 1 0,900 
nautical  mile  altitude).  The  satellites  will  be  uniformly  spaced  (1 20°  apart)with  three  satellites 
in  each  of  six  xbital  planes,  each  xbital  plane  inclined  at  55°  to  the  equetx.  Each  xbital  plane 


will  be  separated  from  each  other  in  longitude  by  60°.  Relative  phasing  of  the  satellites  from  me 
orbital  plane  to  the  next  is  40°,  which  simply  means  that  when  an  ascending  satellite  in  one  plane 
is  crossing  the  equatorial  plane,  an  ascending  satellite  in  the  adjacent  plane  to  the  east  is  40° 
above  the  equatorial  plane  in  its  own  orbital  plane  (15:E9.3. 1).  Each  satellite  transmits 
navigation  signals  at  frequencies  of  1 575.4  MHz  and  1 227.6  MHz,  which  contain  navigation  data 
such  as  satellite  ephemeris  and  satellite  clock  bias  information  ( 2:227).  Use  of  the  two  signals 
permits  the  user's  equipment  to  compensate  for  the  ionospheric  group  delay  or  electromagnetic 
disturbances  in  the  atmosphere  which  may  alter  the  affected  signals  (19:  21.1.1).  Each  satellite 
has  a  mean  mission  duration  of  6.2  years  and  a  design  life  of  7.5  years  (15:  E9.3.7),  and  the 
present  design  provides  for  the  placement  of  three  additional  satellites  within  the  constellation  to 
act  as  active,  in-orbit,  replacement  spares. 

Control  Segment.  The  Ground  Control  Segment  monitors  the  broadcast  satellite  signals 
and  uplinks  corrections  to  ensure  predefined  accuracies.  The  operational  control  segment  will 
consist  of  five  monitor  stations,  a  master  control  station,  and  three  uplink  antennas.  The  widely 
separated  monitor  sets,  positioned  worldwide,  will  allow  simultaneous  tracking  of  the  full  satellite 
constellation  and  will  relay  orbital  and  clock  information  to  the  master  control  station  (2:226). 
The  ranging  data  accumulated  by  the  monitor  stations  will  be  processed  by  the  Navstar  Operations 
Center  (master  control  station)  for  use  in  satellite  determination  and  systematic  error 
elimination  ( 1 9:  2 1 . 1 . 1 ).  The  master  control  station  then  forms  corrections  that  are  uploaded  to 
the  satellites  by  the  uplink  antennas  ( 2:227). 

User  Segment.  The  User  Segment  selects  the  four ,  best  positioned  satellites  from  those 
visible  and,  using  the  navigation  signals  passively  received  from  each  of  these  satellites,  the 
user's  receiver  measures  four  independent  pseudoranges  and  pseudorange  rates  to  the  satellites. 
The  receiver/processor  then  converts  these  signals  to  three-dimensional  position,  velocity,  and 
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Concept  of  Operation 

With  1 6  satellites,  if  the  user  is  able  to  receive  the  navigation  signal  from  satellites  above  a 
five  degree  elevation  angle,  there  are  always  four  or  more  satellites  available  for  use.  ( The  limit 
on  elevation  angle  is  due  to  terrain  masking,  line  -of-sight  obstructions,  etc.)  GPS  satellites 
continuously  transmit  navigation  signals  at  two  L-band  carrier  frequencies:  L I ,  at  1575.42  MHz 
and  L2,  at  1 227.6  MHz.  L 1  is  modulated  by  a  Coarse/Aquisition  (or  C/A)  code  and  a  Precision  (or 
P)  code,  whereas  L2  contains  either  the  P  or  C/A  code  but  not  both;  the  P-code  is  normally  used. 
In  addition,  superimposed  on  both  frequencies  is  system  data  at  50  bits  per  second  that  provide 
handover  information  from  the  C/A  to  P-code,  satellite  orbital  characteristics,  satellite  health 
satatus,  and  satellite  clock  errors  (2:227).  Comparison  of  the  delay  between  the  two  signals 
allows  for  proper  computation  of  errors  due  to  ionospheric  propagation  or  electromagnetic 
disturbances  in  the  atmosphere  ( 1 9:2 1 . 1.1 ). 

Navigation  using  GPS  is  accomplished  by  passive  triangulatian.  To  obtain  a  navigation 
solution,  the  user  measures  "pseudorange"  to  each  of  the  four  satellites  by  timing  the  received 
pseudorandom  code  epoch  with  respect  to  his  local  estimate  of  time.  The  term  “pseudorange"  is 
used  since  the  timing  measurement  of  true  range  to  the  satellite  contains  an  error  in  the  form  of  a 
yet  undetermined  clock  bias.  The  position  of  the  four  satellites  is  computed  using  the  received 
ephemeris  data.  Using  this  data  and  the  pseudorange  measurements,  the  user  then  solves  four 
equations  in  four  unknowns  to  determine  his  three-dimensional  position  and  precise  time. 
Three-dimensional  velocity  is  determined  by  measuring  the  doppler  shifts  on  the  received  carrier 
frequency  (2:227). 


III.  Navigation  Equations 


The  Navigation  Solution 

In  order  to  develop  the  navigation  equations  used  in  solving  for  the  user's  three-dimensional 
position  and  time,  we  will  utilize  the  earth-centered,  right-handed,  Cartesian  coordinate  system 
as  illustrated  in  Figure  A.  At  time  zero,  the  X-axis  passes  through  the  intersection  of  the 
equatorial  plane  and  prime  meridian,  the  Z-axis  passes  through  the  North  Pole,  and  the  Y-axis 
completes  the  right-handed  orthogonal  system.  Because  of  the  earth's  rotation,  the  x  and  y 
coordinates  are  constantly  changing  in  longitude  with  time. 


The  basic  nonlinear  equations  using  four  satellites  are 


Kx-Xj)2  +  (y-yj)2  ♦  (z-2j)2  ft  ♦  T  =  ( I) 

where 

Xj ,  yj ,  Zj  =  position  of  the  ith  satellite  ( known) 

x,y,2  =  user  position  (unknown) 

T  -  clock  bias  ( unknown) 

R{  =  pseudo  range  measurement  to  the  ith  satellite 


The  quantities  Rj.R2.R3  and  R4  are  “pseudoranges“  in  that  they  are  the  sum  of  the  actual 

range  displacements  plus  the  offset  due  to  user  time  error.  Although  these  equations  can  be  solved 
directly,  user  equipment  employs  a  much  simpler  version  of  these  equations  which  can  be  derived 
by  linear i2ation  methods  as  follows  ( 1 4:97): 

Let 


xn,Vn(2n^n 
Ax,  Ay,  A2,  AT 

Rni 


nominal  (8  priori  best  estimate)  values  of  x,y,2,  and  T 

corrections  to  nominal  values 

nominal  pseudorange  measurement  to  ith  satellite 


AR|  =  residual  ( difference)  between  actual  and  nominal 
range  measurements 


Thus  we  obtain  the  following  incremental  relationships: 
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X 


=  xn  +  Ax 

y  -  yn  ♦  Ay 
z  =  zn+Az 
T  =  Tn  +  AT 

*  ^ni * 
and 

R*  - 1< v«.)2  *(y»-y(>2  *<W2)*  *  t. 

Substituting  these  incremental  relationships  into  ( I )  we  obtain 

[(x„  ♦  Ax-  Xj)2  ♦  (yn  ♦  Ay-  y,)2  ♦  (z„  ♦  A  2-  Zj  )2 

=  «ni  +  Aft,  -  Tn  -  AT,  (2) 

Working  with  the  left  hand  side  of  the  equation  and  expanding  terms  we  obtain 

[(xn  ♦  Ax-  x,)2  ♦  (yn  ♦  Ay-  yj)2  *(zn  +  A  z-  z{  )2  )* 

-  f(xn  -  Xj)2  ♦  2  Ax(xn  -x^  ♦  (Ax)2  ♦  (yn  -  y4)2  +  2  Ay(yn  -  y}) 

+  ( Az)2  ♦  (2„  -  z,)2  ♦  2  A 2(2,,  -  z,)2  +  (Az)2]*4  (3) 

Rearranging  terms  and  eliminating  second  order  terms,  this  expression  can  be  written 
simply  as 

(a+ 2b)*4 
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a  =  [(xn-xi)2  +  (yn-yj)2  +  (2n-zj)2] 

b  -  [  Ax( -Xj)  ♦  Ay(yn  -yj)  +  Az(z„  -Zj)] 
Expanding  this  form  using  the  binomial  series  expansion  we  obtain 


(a*  2b)*4  =  (a)14  [J  +  2b/a]!4 

=  (a)l,4[1  ♦  b/a  +  Higher  order  terms] 

By  noting  that  all  higher  order  terms  contain  second  order  and  higher  terms  of  Ax,  Ay,  or 
Az,  we  may  ignore  them  and  our  equation  reduces  to 

(a)54!!  +  b/a]  =  (a)54  ♦  (a^b/a] 

=  (a)54  ♦  bAa)*4 

=  Rni  +  ARj  -  Tn  -  AT  (4) 

Substituting  our  incremental  relationship  for  Rn,  into  this  equation  and  simplifying 

(a)54  ♦  b/ta)14  =  Ka)54*  Tn]  ♦  AR(  -  Tn  -  AT  (5) 


b/ta)*4  =  ARj  -  AT 


R-;  =  (a)*4  ♦  T„  (a)14  =  R»:  -  T» 


Substituting  this  expression  and  our  expression  for  b  into  equation  (6)  we  obtain  the 
linearized  equations  (i=  1,2  3,4)  that  relate  pseudorange  measurements  to  the  desired  user 
navigation  information  as  well  as  the  user's  clock  bias: 

((xn-Xj)/(Rnj  -  Tn)l  Ax  -  l(yn  -Yj)/(Rnj  -  Tn)]  Ay 

+  l(zn-zj)/(RnrTn)l  Az  +  AT 

=  AR,  (7) 

The  quantities  on  the  right-hand  side  of  equation  (7)  are  known;  they  are  simply  the 
differences  between  the  actual  measured  pseudoranges  and  the  predicted  measurements  which  are 
supplied  by  the  user's  computer  based  on  knowledge  of  the  satellite  position  and  current  estimate 

of  the  user's  position  and  clock  bias.  The  quantities  to  be  computed  (  Ax,  Ay,  Az,  and  AT  )  are 
the  corrections  that  the  user  will  make  to  his  current  estimate  of  his  position  and  time  bias.  The 
coefficients  of  these  quantities  on  the  left-hand  side  represent  the  direction  cosines  of  the  line  of 
sight  vector  from  the  user  to  the  satellite  as  projected  along  the  x,  y,  and  z  coordinate  axis 
system  ( 14:98). 

The  four  linearized  equations  represented  by  ( 7)  can  be  expressed  in  matrix  notation  as 


B,2 

*13  1 

A  X 

tR, 

*21 

15*1 

8  22 

1532 

15  23  1 

1533  1 

X 

A  y 

A  2 

= 

A  R._, 

A 

J34, 

J'i42 

1543  1 

AT 

iR4 

_ 

. 

. 

. 

1 


where  Gy  is  the  direction  cosine  of  the  angle  between  the  range  to  the  ith  satellite  and  the  jth 
coordinate  ( M:  98). 

To  express  this  equation  more  compactly  we  will  let 

r  =  the  four  element  pseudorange  measurement  difference  vector 
x  =  user  position  and  time  correction  vector 
B  =  the  A  x  A  solution  matrix 


b11  Bt2  b13  1 

G  n «  Si  >-,r,  \j  tyj  1 

J^31  B  32  #33  ' 

J'i41  G42  £43  1 


x  =  [  a  x  AIJ  a  z  at] 

r  =  l  aR(  aR2  aR3  aR4  1  T 

Thus  our  equation  becomes  simply 

Bx  =  r  or  x  =  B-’  r  (8) 

which  compactly  expresses  the  relationship  between  pseudorange  measurements  and  user  position 
and  clock  bias. 

To  understand  how  the  geometry  of  the  satellites  at  a  point  in  time  can  result  in  a  system 
outage,  we  need  only  examine  the  solution  matrix.  If  the  ends  of  the  unit  vectors  from  the  user  to 
the  four  satellites  selected  are  in  a  common  plane,  the  direction  cosines  of  the  four  unit  vectors 


along  a  direction  perpendicular  to  this  plane  are  all  equal.  When  this  occurs,  the  determinant  of 
the  4  x  A  solution  matrix  becomes  zero  (solution  matrix  becomes  singular)  and  no  solution  is 


possible  from  the  four  equations.  Consequently,  the  navigation  equations  “blow  up"  and  what  is 
known  as  a  "system  outage"  occurs  due  to  the  poor  geometry.  The  situation  where  this  occurs  is 
very  close  to  where  the  four  satellites  are  in  a  common  plane  in  space  as  shown  in  Figure  5  below. 


Figure  5.  i .  -ustration  of  Navstar  Satellite  Geometry  ( 1 1 .433) 


User  Position  Accuracy 

In  order  to  determine  the  accuracy  available  from  the  four  satellites  selected  as  a  function  of 
their  geometry,  we  must  calculate  the  dilution  of  precision  (DOP)  values  available  from  the  four 
satellites  selected.  Since  the  overall  position  accuracy  is  a  product  of  this  value  and  other  system 
errors,  small  DOP  values  are  highly  desirable  in  order  to  yield  accuracies  on  the  order  of  those 
previously  shown  in  Table  I.  We  have  already  seen  how  a  poor  geometry  of  the  four  satellites 
selected  results  In  the  "blowing  up"  of  the  navigation  equations;  the  corresponding  effect  on  the 
DOP  values  is  to  cause  them  to  become  Infinitely  large  (resulting  in  a  system  "outage")  which 
lasts  until  the  geometry  of  the  four  satellites  improves  or  an  additional  satellite  becomes  available 
providing  a  new  combination  with  more  favorable  geometry. 

Since  the  navigation  equation  (8)  we  derived  in  the  last  section  is  a  linear  relationship,  it 
can  also  be  used  to  express  the  relationship  between  errors  in  pseudorange  measurement  and  the 
errors  in  user  position  and  clock  bias.  Mathematically,  the  relationship  can  be  expressed 


Cx 

= 

b->  er 

where 

er 

■ 

pseudorange  measurement  errors 

Cx 

■ 

corresponding  error  in  user  position  and  clock  bias 

If  we  let 

Cr 

■ 

covariance  matrix  of  errors  in  pseudorange  measurements 

and 

Cv 

> 

covariance  matrix  of  the  resulting  errors  in  the  three 
components  of  user  position  and  clock  bias 

then  the  matrix  relationship  between  the  two  covariance  matrices  becomes 


where  6  is  the  same  4x4  matrix  of  coefficients  of  the  unknowns  that  we  derived  previously  and  is 
a  function  only  of  the  direction  cosines  of  the  LOS  unit  vectors  from  the  user  to  the  four  nevsats 
and  the  user's  clock  bias  ( 1 7: 1 0).  Thus,  the  error  relationships  are  functions  only  of  satellite 
geometry,  which  leads  to  the  concept  of  geometric  dilution  of  precision  as  a  measure  of  how 
satellite  geometry  degrades  accuracy  ( 1 4:99). 

In  order  to  provide  a  method  of  quantitatively  determining  whether  a  given  satellite 
geometry  is  good  or  bad,  we  will  assume  that  the  geometric  effects  on  the  errors  in  pseudorange 
measurements  are  very  small  (a  reasonable  assumption).  Thus  ,  a  good  approximation  of  the 
geometric  performance  can  be  obtained  by  letting  each  individual  pseudorange  measurement  have  a 
one  sigma  error  of  unity;  thus,  the  covariance  matrix  for  the  errors  in  pseudorange 
measurements  becomes  a  4  x  4  identity  matrix  ( 1 4:99).  This  means  that  the  ranging  errors  from 
each  satellite  are  assumed  to  be  unity,  completely  random,  and  that  the  correlation  of  ranging 
errors  between  satellites  is  zero.  To  a  good  approximation, 

Cy  =  (BtB)-', 

and  assuming  sufficient  signal  strength,  this  covariance  matrix  depends  only  on  the  direction 
and  is  in  no  way  dependent  on  the  distances  between  the  user  and  each  satellite  ( 1 7: 1 0).  Since 
the  diagonal  elements  of  this  covariance  matrix  are  actually  the  variances  of  user  position  and 
time,  various  OOP  values  can  be  obtained  from  looking  at  appropriate  elements  of  this  matrix. 
These  values  can  then  be  used  to  compare  position  accuracies  (from  a  geometry  standpoint)  of 
different  orbital  configurations  as  well  as  to  measure  the  overall  effect  of  geometry  to  errors  in 
the  user  s  position  and  time  bias. 


Dilution  of  Precision  Values 


If  we  express  the  diagonal  elements  of  the  covariance  matrix  of  errors  in  user  position  and 
clock  bias  as  Yx,  Yy,  V2,  and  VT,  then  the  four  dimensional  geometrical  dilution  of 
precision  ( GDOP )  is  obtained  by  taking  the  square  root  of  the  trace  of  the  matrix: 

GDOP  *  [  Vx  ♦  Yy  +  V2  ♦  VT)* 

This  factor  includes  all  four  unknowns  (three  dimensions  of  position  and  time)  and  is  the 
conventional  measure  of  overall  geometric  performance. 

A  more  frequently  used  measure  of  geometric  performance  is  the  three  dimensional 
position  dilution  of  precision  (PDOP)  which  relates  only  to  the  three  components  of 
position  error.  PDOP  is  also  invariant  with  the  coordinate  system  and  is  used  because  the  most 
important  consideration  in  any  navigation  system  is  position  accuracy;  knowing  time  is  a 
secondary  byproduct  ( M:100).  It  is  defined  as  the  square  root  of  the  sum  of  the  squares  of  the 
three  components  of  position  error ,  or  mathematically: 

PDOP  =  (Yx  +  Yy  ♦  Y2]* 

Several  other  alternative  OOP  values  are  also  occasionally  used  in  evaluating  satellite 
constellations  and  relate  only  some  of  the  variances  of  user  position  and  time.  These  include  the 
horizontal  dilution  of  precision  (HDOP),  the  altitude  dilution  of  precision  (VDOP), 
the  time  dilution  of  precision  (TDOP),  and  the  larger  component  of  horizontal 
position  error  ( MDOP ).  HDOP ,  VDOP ,  TDOP ,  and  MDOP  ere  the  DOP  factors  that  apply  for  the 
horizontal  position  error,  the  altitude  error,  the  error  in  user  clock  bias  (multiplied  by  the 
speed  of  light),  and  the  larger  component  of  the  horizontal  position  error ,  respectively.  They  ire 
defined  mathematically  as  follows  ( 1 7: 1 0): 


HOOP 


*  [V«  *  V14 

YDOP  =  [Vz]^ 

TDOP  =  [Vt-,14 

MDOP  =  Max[(Vx)^,(Vy)^] 

Satellite  Selection 

In  order  to  obtain  the  most  accurate  user  position,  it  would  be  highly  desirable  to  utilize 
those  four  satellites  with  the  most  favorable  geometry  (lowest  DOP  values)  with  respect  to  the 
user  at  any  instant  of  time.  This  presents  no  problem ,  should  there  be  only  four  visible  satellites 
to  choose  from .  as  all  four  must  be  used  to  determine  the  user's  three  dimensional  position.  The 
majority  of  the  time,  however ,  there  will  be  six  or  more  satellites  In  view  by  an  earth  based  user 
and  even  more  by  a  low  altitude  satellite  user,  and  the  computet ionel  time  required  to  compute 
POOP  values  for  all  the  possible  combinations  of  satellites  is  excessive. 

The  results  of  many  computer  runs  and  analytical  studies  have  demonstrated  an  almost  total 
correlation  between  POOP  and  the  volume  of  a  tetrahedron  formed  by  lines  connecting  the  tips  of 
the  four  unit  vectors  from  the  user  toward  the  four  navsats  (17  11).  Usually,  the  larger  the 
volume  of  this  tetrahedron,  the  smaller  the  corresponding  PDOP  value  will  be  for  this  same  set  of 
satellites.  Since  the  computational  time  for  computing  the  volume  of  a  tetrahedron  for  each 
different  combination  of  satellites  is  much  less  than  the  time  required  to  calculate  a  POOP  value 
(which  requires  a  matrix  inversion),  the  computer  program  used  in  this  work  as  well  as  many 
other  similar  studies  in  the  pest  is  designed  to  first  compute  the  volumes  of  the  tetrahedrons 
associated  with  each  different  combination  of  four  satellites,  identify  the  "best  four"  which  yield 
the  largest  tetrahedron  volumes,  and  then  use  that  combination  of  satellites  to  compute  the  DOP 


20 


values  ( 1 7: 1 2).  Other  methods  of  selecting  the  optimum  four  satellites  may  possibly  be  more 
efficient  and  have  been  the  subject  of  many  analytical  studies,  but  they  will  not  be  addressed  in 
this  paper.  The  problem  of  satellite  selection  is  not  only  a  problem  to  those  conducting  global 
analyses  of  orbital  constellations  for  evaluation,  but  it  is  just  as  significant  a  problem  to  the 
designers  of  User  Set  Equipment,  as  the  equipment  must  be  designed  to  operate  quickly  tfn  a 
dynamic  environment  and  large  computer  resources  are  not  available. 
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IV.  Baseline 


Analysis 


Global  Distribution  Bun 

Introduction.  The  reference  orbit  parameters  for  the  OPS  baseline  constellation  used  in 
this  analysis  are  given  in  Table  II.  The  computer  program  used  in  conducting  the  analysis  was  a 
modified  version  of  a  program  developed  by  the  Rand  Corporation  of  Santa  Monica,  California, 
on  the  geometric  performance  of  pseudoranging  navigation  satellite  systems  ( 1 7).  This  program , 
which  itself  is  a  modification  of  an  earlier  program,  developed  by  the  Aerospace  Corporation 
for  earth-based  users,  was  selected  for  the  evaluation  over  the  computer  program  currently  being 
used  by  the  Aerospace  Corporation,  EGAD  (Efficient  GPS  Availability  Determination)  Program 
(5:2),  for  several  reasons.  First,  the  Rand  Corporation  program  had  been  modified  to 
accommodate  users  in  any  earth  orbit,  while  the  EGAD  program  did  not  have  any  provision  for  the 
space-based  user.  The  EGAD  program,  although  a  much  longer  program  with  many  more 
capabilities,  restricted  the  shape  and  size  of  the  orbits  of  the  navigation  satellites;  the  program 
modified  by  the  Rand  Corporation  allows  the  user  to  select  any  size  or  shape  of  the  orbits  of  either 
the  navsats  or  the  user  satellite.  Finally,  the  Rand  Corporation  program  Includes  a  feature  which 
allows  for  the  variation  of  the  navsat  antenna  beamwidth  and  determines  the  effect  of  this  variation 
on  navigational  accuracy  for  satellite  users;  the  EGAD  Program  has  no  such  capability. 

A  modification  was  made  to  the  Rand  program  to  allow  for  the  determination  of  the  effect  of 
satellite  losses  on  the  baseline  or  any  other  alternative  constellation,  tn  addition,  a  few 
corrections  were  found  necessary  in  order  for  the  program  to  give  reliable  statistics  on  DOP's 
when  fewer  than  four  satellites  were  visible.  A  minor  modification  of  the  program  output  was  also 
made  to  provide  the  number  of  satellites  visible  to  the  space-based  user  when  fewer  than  four 
satellites  were  available  to  the  user  satellite.  A  description  of  the  program  and  its  modifications 


is  provided  in  Appendix  B,  and  the  complete  computer  program  listing  is  given  in  Appendix  C. 


TABLE  II 

Baseline  Constellation  Orbital  Parameters  ( 13:  D2.3.3) 


Satellite  Orbit  Longitude  of  the  Right  Ascension  of 

Number  Plane  Ascending  Node  f  Deal  the  Ascending  Node  [Deo8! 


1 

1 

000 

030 

2 

1 

060 

030 

3 

1 

120 

030 

4 

2 

080 

090 

5 

2 

140 

090 

6 

2 

020 

090 

7 

3 

160 

150 

8 

3 

040 

150 

9 

3 

100 

150 

10 

4 

060 

210 

11 

4 

120 

210 

12 

4 

000 

210 

13 

5 

140 

270 

14 

5 

020 

270 

15 

5 

080 

270 

16 

6 

040 

330 

17 

6 

100 

330 

18 

6 

160 

330 

•  Referenced  to  astronomical  coordinates  of  1 950.0  as  of  1  July  1 985, 
0  hr  0  min  GMT  and  regressing  at  -0.04009  deg/d ay. 


Selection  Parameters  For  all  global  distribution  runs,  a  uniform  distribution  of 
users  is  approximated  by  the  OOP’s  of  users  at  a  given  latitude  by  the  cosine  of  that  latitude. 
Whenever  the  constellation  selected  was  a  symmetrical  arrangement  of  satellites,  only  the  OOP’s 


for  the  northern  hemisphere  were  calculated  as  this  permitted  a  larger  number  of  sample  points 
to  be  selected  for  the  analysis.  Since  the  arrangement  is  symmetrical,  an  analysis  (on  a  global 
basis  only)  of  the  same  number  of  sample  points  uniformly  distributed  in  the  southern 
hemisphere  would  yield  identical  results.  Consequently ,  the  results  shown  for  the  analysis  of  the 
northern  hemisphere  are  statistically  representative  of  the  entire  globe. 

The  proper  latitude  and  longitude  step  sizes  were  determined  from  the  analysis  of  several 
data  runs  using  alternative  parameters.  The  largest  step  sizes  that  would  provide  statistically 
representative  values  were  selected  to  provide  a  good  balance  between  accuracy  and  computer 
computation  costs,  which  are  enormous  for  global  distribution  runs.  The  time  increments  used 
were  selected  in  a  similar  manner. 

For  symmetrical  constellations,  the  time  interval  selected  for  each  run  was  four  hours, 
since  after  this  length  of  time,  each  satellite  within  thr;  ^stellation  would  have  moved  to  the 
original  position  of  the  satellite  adjacent  (within  the  same  orbit)  to  it  at  the  start  of  the  run. 
( Each  satellite  is  separated  from  each  other  within  the  orbital  planes  by  four  hours.)  Sines  this 
analysis  was  concerned  only  with  the  geometric  performance  of  the  baseline  constellation  for  a 
global  basis,  rather  than  for  specific  user  locations,  this  four  hour  time  period  provides 
statistically  good  data  for  any  24  hour  period,  due  to  the  repetiveness  of  the  satellite  motion. 

For  the  non-symmetrical  constellations,  the  time  interval  analyzed  was  increased  to  six 
hours,  or  half  the  1 2  hour  orbital  periods  of  the  GPS  satellites.  This  is  the  minimum  time  period 
that  could  be  chosen  to  provide  statistically  representative  data  (on  a  global  basis)  for  any  24 
hour  period.  As  a  result  of  tho  longer  time  interval  required  for  each  run,  the  time  increments 
were  doubled  in  order  to  remain  within  computer  constraints  and  yet  still  provide  reasonably 
accurate  data.  Since  the  analysis  in  these  cases  focused  on  the  comparative  results  between 
constellations  analyzed,  rather  than  exact  OOP  values,  this  step  size  was  found  to  be  adequate. 


Baseline  Constellation  Orbital  Elements 


ORBITAL  ELEMENTS 


ECC 

ARGP 

(Deg) 

RASC 

(Deg) 

INC 

(Deg) 

ANOM 

(Deg) 

1 

0.00 

0.00 

30.00 

55.00 

0.00 

2 

0.00 

0.00 

30.00 

55.00 

120.00 

3 

0.00 

0.00 

30.00 

55.00 

240.00 

4 

0.00 

0.00 

90.00 

55.00 

40.00 

5 

0.00 

0.00 

90.00 

55.00 

160.00 

6 

0.00 

0.00 

90.00 

55.00 

280.00 

7 

0.00 

0.00 

150.00 

55.00 

80.00 

8 

0.00 

0.00 

150.00 

55.00 

200.00 

9 

0.00 

0.00 

150.00 

55.00 

320.00 

10 

0.00 

0.00 

210.00 

55.00 

120.00 

11 

0.00 

0.00 

210.00 

55.00 

240.00 

12 

0.00 

0.00 

210.00 

55.00 

360.00 

13 

0.00 

0.00 

270.00 

55.00 

160.00 

14 

0.00 

0.00 

270.00 

55.00 

280.00 

15 

0.00 

0.00 

270.00 

55.00 

40.00 

16 

0.00 

0.00 

330.00 

55.00 

200.00 

17 

0.00 

0.00 

330.00 

55.00 

320.00 

18 

0.00 

0.00 

330.00 

55.00 

80.00 

(I 


PARAMETERS  USED  IN  GLOBAL  DISTRIBUTION  CALCULATIONS: 

MASKING  ANGLE  -  5.00  DEGREES 
LATITUDE  STEP  -  10.00  DEGREES 
LONGITUDE  STEP  -  20.00  DEGREES 
TOTAL  TIME  (MIN)  «  240 
TIME  INCREMENT  (MIN)  -  5 
HEMISPHERE  EVALUATED  -  NORTHERN 


Results.  If  we  define  an  "outage"  to  be  situations  where  the  POOP  for  a  given  location  8t 
some  time  exceeds  6.0  ( a  1  -sigma  ranging  error  of  7  meters  would  provide  a  RMS  error  of  6.0  x 
7m,  or  42  meters  in  this  case),  as  is  commonly  chosen  on  analyses  of  this  nature,  then  we  can 
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This  agrees  closely  with  a  similar  analysis  conducted  earlier  by  the  Aerospace 
Corporation  on  the  same  six-plane  constellation  ( 1 5:  E9.3.2) ,  and  verifies  the  accuracy  of  the 
data  If  we  choose  to  alternatively  define  an  outage  as  a  PDOP  greater  than  7.0,  as  some  studies  on 
earlier  constellation  designs  have  assumed,  then  we  can  see  that  PDOP's  greater  than  7.0  occur 
only  0.41  X  of  the  time  for  the  proposed  baseline  constellation.  Although  this  paper  is  concerned 
primarily  with  the  PDOP  values,  statistics  for  all  six  DOP's  are  listed  in  Table  IV,  and  a  complete 
breakdown  of  each  OOP  by  latitude  is  provided  in  the  appendix. 

The  maximum  and  minimum  number  of  satellites  available  to  the  user  at  each  latitude  and 
longitude  are  shown  in  Table  V.  (Note  that  not  all  latitudes  and  longitudes  are  observed.)  The 
probabilities  of  n  or  more  satellites  b8ing  visible  above  a  5  degree  elevation  angle  are  shown  in 
Table  VI.  For  the  proposed  1 8-satellite  baseline  constellation,  at  least  four  satellites  are  always 
available  to  the  earth-based  user,  and  consequently,  the  outages  that  occur  ore  due  solely 
to  the  poor  geometry  of  the  satellites  available  to  the  user.  By  examining  the 
computer  output  located  in  the  appendix  on  each  individual  DOP  for  each  latitude,  we  can  also 
determine  the  location  of  these  outage  areas  during  any  24  hour  period.  Typical  location  of  outages 
and  time  durations  for  the  six-plane  constellation  are  depicted  in  Figure  6, 7, 8,  and  9.  As  one  can 
observe  from  either  these  pictorial  presentations  or  by  the  computer  output  located  in  the 
appendix ,  the  primary  outage  locations  for  the  baseline  constellation  occur  in  pairs  approximately 
centered  on  latitudes  of  35N  and  65S,  with  corresponding  pairs  centered  at  35S  and  65N ,  as  shown 
in  Figure  7. 
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TABLE  V 


Maximum  and  Minimum  Numbers  Seen  at  Each  Latitude  and  Longitude 


MAXIMUM 

LONGITUDE  (DEG) 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

3 

3 

3 

2 

4 

6 

8 

0 

2 

4 

6 

8 

0 

2 

4 

6 

8 

0 

2 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

7 

8 

8 

8 

8 

7 

8 

8 

8 

7 

8 

8 

8 

8 

7 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

6 

7 

8 

7 

6 

7 

8 

8 

7 

6 

7 

8 

7 

6 

7 

8 

8 

7 

6 

7 

8 

7 

6 

6 

8 

8 

7 

6 

7 

8 

7 

6 

6 

8 

8 

7 

7 

7 

7 

6 

7 

7 

7 

7 

7 

7 

7 

7 

6 

7 

7 

7 

7 

7 

8 

8 

8 

7 

8 

8 

8 

8 

7 

8 

8 

8 

7 

8 

8 

8 

8 

7 

8 

7 

8 

8 

8 

8 

7 

8 

7 

8 

7 

8 

8 

8 

8 

7 

8 

MINIMUM 


LONGITUDE  (DEG) 


TABLE  VI 


LATITUOI 
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System  Degradation  Due  to  satellite  Losses 


Introduction.  The  effect  of  satellite  losses  from  the  1 8-satellite  baseline  constellation 
was  also  analyzed  on  a  global  basis,  as  this  was  one  of  several  factors  considered  in  selecting  the 
optimum  configuration  of  the  1 8  satellites.  For  this  part  of  the  analysis,  it  was  assumed  that  no 
active,  in-orbit  spares  were  available.  The  effect  upon  geometric  performance  from  losses  of  one 
to  three  satellites  was  evaluated  was  evaluated  for  several  representative  cases,  with  an  emphasis 
on  determining  the  best  and  worse  cases  possible,  instead  of  on  the  expected  performance  due 
to  random  satellite  losses,  as  most  earlier  studies  have  dealt  with.  "Best  case"  information 
might  be  particularly  useful  to  those  involved  in  the  initial  buildup  of  the  constellation  to  its  full 
operational  1 8-satellite  configuration  as  a  means  of  determining  both  the  order  and  placement  of 
each  successive  satellite  deployment  Likewise,  "worst  case"  scenarios  might  prove  to  be  useful  to 
those  military  planners  seeking  defenses  against  anti-satellite  weapons  that  potential  enemies 
might  design  and  deploy.  Since  computer  restraints  did  not  permit  all  possible  combinations  of 
satellite  losses  to  be  analyzed  (for  the  three  satellite  case  there  are  8 1 6  possible  combinations), 
the  satellite  losses  analyzed  were  carefully  selected  in  order  to  determine  the  most  probable 
best  and  worst  cases.  In  addition,  marry  of  the  combinations,  due  to  the  symmetric  arrangement  of 
the  satellites  and  their  repetitive  nature,  provided  identical  degradation  of  performance  (from  a 
global  standpoint)  and  could  be  excluded.  For  example,  If  only  one  satellite  is  lost  (or  destroyed), 
any  satet/fte  selected  will  degrade  the  geometric  performance  in  the  same  way  as  any  other 
one  satellite  lost.  (Only  the  location  of  the  outages  will  differ,  which  could  be  of  interest  to 
military  planners).  Similarly,  a  loss  of  a  pair  of  satellites  (1)  and  (4)  degrade  system 
performance  the  same  amount  (on  a  global  basis)  as  a  loss  of  satellites  ( 2)  and  (5),  (3)  and  (6), 
( 5)  and  ( 8),  and  many  other  symmetrically  equivalent  pairs.  ( One  need  only  realize  that  in  four 
hours  satellites  ( 1 )  and  ( 4)  move  to  the  same  inertial  positions  as  ( 2)  and  ( 5)  were  originally 


located  in  order  to  understand  this). 


Results.  By  examining  the  tabular  data  on  some  representative  runs  in  Table  VII ,  one  can 
easily  see  that  which  satellites  are  deleted  from  the  constellation  greatly  affect  the  amount  of 
degradation  to  system  performance. 


TABLE  VII 

System  Performance  Degradation  Due  to  Satellite  Losses 


Satellites 

%  POOP 

%  POOP 

Deleted 

>  6-0 

>LG 

None 

.0055 

.0041 

6 

.0230 

.0198 

1 

.0232 

.0201 

1.5* 

.0478 

.0407 

1.12 

.0488 

.0458 

1.2 

.0507 

.0434 

1.7 

.0515 

.0450 

1.8 

.0720 

.0634 

1,6 

.0749 

.0695 

1.4 

.0779 

.0697 

1.11 

.0782 

.0716 

1.9 

.0787 

.0709 

1,10** 

.0800 

.0725 

4.5.6* 

.0780 

.0679 

1,7.13 

.1017 

.0900 

1.7.8 

.1042 

.0920 

1,4.5 

.1057 

.0928 

1,2.4 

.1349 

.1245 

1,4.7 

.1380 

.1258 

1,9,15 

.1387 

.1289 

1,10.11 

.1435 

.1325 

1,8,15** 

.1584 

.1454 

•Best  Case 
**  Worst  Case 


As  is  evident  from  this  table,  the  percentage  of  time  a  system  outage  occurs  in  the  two 
satellite  loss  “worst  case"  scenario  is  actually  greater  than  the  percentage  of  time  that  an  outage 
occurs  in  the  three  satellite  loss  “best  case"  (8.00*  versus  7.80*  ).  The  rangB  of  values 
(expressed  as  percentages  of  POOP  >  6)  is  rather  widely  spread  from  the  best  to  worst  cases, 
varying  from  4.78*  to  8.00*  in  the  two  satellite  C8se  and  from  7.80*  to  15.84*  in  the  three 
satellite  case.  From  examining  this  data,  it  appears  that  a  loss  of  closely  grouped  satellites,  such 
as  three  adjacent  satellites  within  the  same  plane,  provide  the  "best  case"  situation,  while,  as 
might  be  expected,  losses  of  satellites  widely  separated  generally  result  in  the  ‘worst  case" 
scenario.  For  the  three  satellite  ( loss)  cases  analyzed,  this  is  easily  deduced  from  observing  that 
the  best  case  occurs  when  all  satellite  losses  occur  in  the  same  plane,  while  the  worst  case  occurs 
when  the  setellitesdeleted  are  separated  to  the  maximum.  From  this  information,  it  appears  that 
should  satellite  losses  be  sustained,  a  significant  improvement  in  performence  might  be 
obteined  in  the  interim  (until  replecement  spores  could  be  hunched)  by 
rephasing  the  remaining  satellites  to  provide  as  wide  as  possible  average 
separation  between  all  satellites  per  period.  Should  the  three  satellite  spares  already  be 
on-orbit  as  presently  planned,  it  suggests  that  these  spares  would  be  optimally  positioned  in  every 
other  orbital  plane,  and  repositioned  as  necessary  within  each  orbit  when  losses  are  sustained,  to 
obtain  the  widest  possible  distribution  of  remaining  satellites.  Presently,  the  planned  location  for 
placement  of  these  three  additional  in-orbit  spares  is  depicted  in  Figure  10,  and  should  provide 
this  opportunity. 
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Satellite  User 


Introduction.  Although  th8  Global  Positioning  System  was  designed  for  use  by  the  earth 
or  near-earth  based  user,  its  potential  application  for  autonomous  navigation  of  satellites  and 
space  navigation  is  highly  desirable.  Until  recently,  few  studies  have  been  conducted  in  this  area, 
yet  based  on  the  those  studies  that  have  been  made  regarding  this  potential  application,  the  results 
appear  promising.  At  least  one  analysis  has  shown  that  a  three-dimensional  accuracy  on  the  order 
of  100  meters  can  be  obtained  for  autonomous  navigation  of  geosynchronous  satellites 
( 1 3:D2.3. 1 ) ,  and  even  better  accuracies  can  be  achieved  for  space-  based  users  at  lower  altitudes. 

In  analyzing  the  geometric  performance  of  the  baseline  constellation  for  the  satellite  user, 
the  antenna  pattern  design  of  the  GPS  satellites  severely  limits  the  times  and  numbers  of  GPS 
satellites  that  can  be  observed  by  the  satellite  user,  particularly  at  the  higher  altitudes.  In  this 
somewhat  limited  analysis  of  the  geometric  performance  of  the  GPS  satellites  for  the  satellite 
user,  its  potential  for  navigation  of  the  spaceborne  user  was  examined  by  analyzing  OOP’s 
obtainable  for  a  user  in  low-earth  orbit,  high-eerth  orbit,  and  intermediate  elliptical  orbits,  as 
well  as  for  a  high  altitude  ballistic  missile  trajectory.  The  effect  of  changing  the  antenna  beam 
angle  of  the  GPS  satellites  upon  the  geometric  performance  of  the  system  for  these  spaceborne 
users  was  also  analyzed  to  determine  if  significant  improvement  could  be  obtained,  particularly 
for  the  high  altitude  user. 

Antenna  Design.  The  GPS  satellites  utilize  a  helical  array  that  provides  a  conical  shaped 
beam  to  cover  the  earth  uniformly  from  their  altitude  of  10,900  nautical  miles.  This  antenna 
array  also  provides  a  skirt  within  its  main  lobe  that  extends  beyond  the  edge  of  the  earth  to  a 
limited  degree,  and  it  is  this  part  of  the  antenna  pattern  that  the  geosynchronous  satellite  user 
may  utilize  for  navigation.  For  the  L I  frequency,  the  main  beam  antenna  gain  at  the  edge  of  the 
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earth  is  approximately  1 4  dB ,  dropping  gradually  to  about  3  dB  at  about  21.4  degrees  from  the 
satellite  to  earth  center  line,  or  an  1 1  dB  variability  in  antenna  gain  (Figure  1 1 )  over  this  skirt 
region  ( 1 3:D2.3.2). 


ANGLE  FROM  ANTENNA  CENTER  LINE  IN  DEGREES 
Figure  1 1 .  Main  Beam  Antenna  Gain  at  L 1  Frequency  ( 1 3:D2.3.3) 


As  the  user  spacecraft  approaches  the  GPS  satellites,  the  NAYSTAR  satellites  keep  their 
antennas  directed  toward  the  center  of  the  earth,  thus  forcing  the  user  to  obtain  signals  coming 
from  the  NAYSTAR  sidelobes  and/or  back  lobe  ( 1 1:432),  as  shown  in  Figure  12.  Although  at  high 
altitudes,  users  will  have  only  the  edge  of  the  antenna  pattern  available,  this  could  be  adequate  for 
updating  ephemerides  if  an  accurate  on-board  clock  is  available  ( 20:2 1 .2. 1 ). 


Low  Altitude  Earth  Orbits.  A  150  nautical  mile  altitude  circular  parking  orbit 
(Case  l )  was  selected  as  a  typical  low  altitude  earth  orbit  for  the  evaluation.  The  orbit  chosen  had 
an  inclination  of  28  degrees  and  a  90  minute  period.  DOP  values  were  computed  at  increments  of  5 
minutes  over  a  total  time  of  12  hours.  For  the  designed  antenna  half-angle  beemwidth  of  21.4 
degrees,  the  POOP  values  obtained  were  extremely  good,  as  were  expected,  and  ranged  from  a  best 
value  of  1.61  to  a  worst  value  of  2.155.  Since  the  PDOP  values  using  the  designed  antenna 
beemwidth  were  significantly  better  than  even  for  the  typical  earth- based  users,  evaluation  at 
other  antenna  beamwidths  was  not  necessary. 

Intermediate  Altitude  Earth  Orbits.  For  this  evaluation,  a  highly  elliptical 
Hohman  transfer  orbit  (Case  2)  from  a  low  altitude  parking  orbit  to  geosynchronous  orbit  was 
chosen  as  representative  of  an  intermediate  altitude  earth  orbit.  The  eccentricity  of  the  selected 
orbit  was  0.73,  and  it  had  a  0  degree  inclination.  The  orbital  period  wss  1 0.56  hours.  DOP  values 
were  computed  at  5  minute  intervals  over  a  12  hour  period.  When  analyzed  with  the  designed 
antenna  half-angle  beemwidth  of  21 .4  deyees,  four  or  more  satellites  were  visible  only  1 7.9*  of 
the  time  ,  three  or  more  satellites  for  24. 1  *  of  the  time,  and  two  or  more  satellites  for  5 1  *  of 
the  time.  No  satellites  were  visible  to  the  user  in  this  orbit  29*  of  the  time. 

The  same  satellite  user  orbit  was  also  analyzed  using  an  antenna  half-angle  beemwidth  of  45 
degrees,  with  significantly  better  results.  With  this  modification  to  antenna  design,  four  or  more 
satellites  were  available  to  the  user  at  all  times,  with  PDOP  values  ranging  from  1 .539  to 
38.663.  PDOP  values  less  than  6.0  were  achieved  35.2*  of  the  time.  The  results  are 


summarized  in  Table  VIII. 


Table  VI 1 1 


Baseline  Constellation  Evaluation  -  Satellite  User 


%  of  T  ime  n  or  More  Satel  1  ites 


User 

Antenna 

Visible 

Satellite 

BWIDTH 

1 

2 

5 

4 

POOP  Range 

Case  1 

21.4 

100.0 

100.0 

100.0 

100.0 

1.61  -2.155 

Case  2 

21.4 

71.0 

51.0 

24.1 

17.9 

1.544- 

45.0 

100.0 

100.0 

100.0 

100.0 

1.539-38.663 

Case  3 

21.4 

100.0 

100.0 

100.0 

100.0 

1.537-  1.818 

Cose  4 

21.4 

62.8 

29.5 

12.4 

7.8 

1.505- 

45.0 

100.0 

100.0 

100.0 

100.0 

1.523-127.57 

90.0 

100.0 

100.0 

100.0 

100.0 

1.504-9.361 

Case  5 

21.4 

55.2 

22.1 

1.4 

0.0 

45.0 

100.0 

100.0 

100.0 

100.0 

14.55-93.023 

90.0 

100.0 

100.0 

100.0 

100.0 

6.199-7.928 

180.0 

100.0 

100.0 

100.0 

100.0 

5.560  -  7.505 

Low  Altitude  1CBM  Trajectory.  A  typical  ICBM  profile  (Case  3)  was  chosen  for 
this  portion  of  the  evaluation.  The  elliptical  trajectory  selected  had  an  eccentricity  of  .55  and  a 
period  of  about  1 .2  hours.  Time  of  flight  from  launch  to  Impact  was  approximately  40  minutes, 
and  the  maximum  altitude  attained  was  slightly  less  than  1 400  nautical  miles.  OOP  values  were 
computed  at  one  minute  intervals  for  the  duration  of  the  flight  As  in  the  low  altitude  earth  orbit 
(Case  1 ),  exceptionally  good  POOP  values  were  achieved  with  the  designed  antenna  beamwidth 
throughout  the  flltfit,  ranging  from  1.537  to  1.818,  as  shown  in  Table  VIII.  Four  or  more 


satellites  were  always  visible,  and  the  POOP  values  improved  with  increasing  altitude  all  the  way 
to  approximately  1 400  nautical  miles. 

High  Altitude  Ballistic  Missile  Trajectory.  A  highly  elliptical  ballistic 
trajectory  (Case  4)  with  an  eccentricity  of  .82  and  a  period  of  1 1.05  hours  was  analyzed. 
Maximum  altitude  attained  during  flight  was  approximately  2 1 ,200  nautical  miles.  DOPs  were 
computed  at  five  minute  intervals  for  the  duration  of  the  flight.  For  the  designed  antenna 
half-angle  beamwidth,  four  or  more  satellites  were  visible  only  7.8$  of  the  time,  three  or  more 
satellites  only  1 2.4  $  of  the  time,  and  two  or  more  visible  only  29.5$  of  the  time.  There  were  no 
satellites  visible  37.2$  of  the  time. 

When  the  same  case  was  analyzed  using  an  antenna  half-angle  beamwidth  of  45  degrees,  four 
or  more  satellites  were  always  visible,  with  PDOP  values  ranging  from  1 .523  to  1 27.569.  POOP 
values  less  than  6.0  were  achieved  20.9$  of  the  time.  As  in  the  previous  case,  PDOP  values 
actually  improved  from  launch  up  to  approximately  1400  to  1800  nautical  miles  in  altitude. 
When  the  antenna  half-angle  beamwidth  was  increased  to  90  degrees,  PDOP  values  improved 
significantly,  ranging  from  1 .504  to  9.36 1 ,  and  PDOP  values  less  than  6.0  were  achieved  50.3$ 
of  the  time,  and  values  less  than  7.0  were  achieved  over  72$  of  the  time. 

High  Altitude  Earth  Orbits.  A  circular,  geosynchronous  earth  orbit  (Case  5)  was 
selected  for  evaluation  and  analyzed  for  several  different  antenna  configurations.  DOPs  were  again 
computed  at  five  minute  intervals  over  a  1 2  hour  period.  With  the  designed  antenna  beamwidth 
there  were  never  four  satellites  visible.  Three  or  more  satellites  were  visible  only  1.4$  of  the 
time,  two  or  more  22. 1  $  of  the  time,  and  one  or  more  visible  only  55.2$  of  the  time.  The 
geosynchronous  user  can  observe  no  satellites  44.8$  of  the  time  wtth  the  designed  antenna 
beamwidth. 

When  the  antenna  half-angle  beamwidth  was  increased  to  45  degrees  (and  higher),  four  or 
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more  satellites  were  always  visible  (Table  VIII),  and  PDOP  values  improved  significantly  with 
increasing  beamwidth  angle.  At  an  antenna  half- angle  beamwidth  of  90  degrees,  PDOP  values  less 
than  7.0  were  achieved  62. 1  %  of  the  time.  When  increased  even  further  to  the  maximum  1 80 
degrees,  PDOP  values  less  than  7.0  were  achieved  87.6*  of  the  time,  ranging  from  values  of  5.56 
to  7.505. 


Figure  13.  Average  Satellite  Visibility  as  a  Function  of  Altitude 

Summary.  For  the  low  altitude  satellite  user ,  the  baseline  OPS  constellation  will  provide 
accuracies  far  exceeding  those  obtained  on  the  earth's  surface,  as  for  these  altitudes  ( less  than 
about  1800  nautical  miles),  PDOP  values  generally  improve  with  increasing  altitude.  More 
satellites  are  generally  visible  to  choose  from,  resulting  in  better  geometry,  ss  shown  in  Figure 
13.  As  the  satellite  user  approaches  higher  altitudes,  however,  the  constraints  due  to  antenna 
beamwidth  cause  a  significant  degradation  of  performance  and  result  in  frequent  outages  caused  by 
the  insufficient  number  of  satellites  available  to  the  user.  Even  in  these  situations,  however, 


navigation  can  be  made  possible  by  means  of  sequential  measurements  from  the  GPS  satellites  as 
they  become  visible  to  the  user ,  8  precise  on-board  clock ,  and  the  knowledge  of  the  known  orbital 
dynamics  of  the  satellite.  Furthermore,  the  improvements  in  accuracy  that  can  be  obtained  by 
increasing  the  antenna  half-angle  beamwidth  can  certainly  not  be  overlooked,  as  such  a 
modification  of  antenna  design  could  significantly  increase  the  position  accuracy  available  to  the 
high  altitude  satellite  user  and  should  be  considered  if  GPS  is  to  be  designed  for  space  navigation  as 
well  as  for  tactical  earth  navigation. 

Another  limitation  on  the  use  of  GPS  for  space  navigation  is  geometrical  in  nature.  Not  only 
is  it  necessary  to  have  four  satellites  in  view  for  three-dimensional  positioning,  but  the  positions 
of  th«se  satellites  should  be  widely  distributed  about  the  user  in  all  directions  to  minimize  PDOP 
values.  Even  with  unlimited  antenna  gain  awl  360  degree  coverage  (180  degree  antenna 
half-angle  beamwidth)  there  is  a  limit  to  the  position  accuracy  that  can  be  achieved  at  higher  and 
higher  altitudes,  ss  the  relative  separation  of  the  GPS  satellites  with  respect  to  the  user  continues 
to  diminish  with  distance.  Consequently,  the  altitude  dilution  of  precision  (VDOP)  becomes  larger 
and  larger,  and  primarily  as  a  result  of  this  growing  component,  the  three  dimensional  position 
dilution  of  precision  (PDOP)  becomes  exceedingly  high.  For  this  reason,  GPS  will  probably  never 
be  practical  for  space  navigation  between  planets  in  the  solar  system  snd  beyond,  but  certainly  is 
feasible  for  satellite  navigation  or  for  the  navigation  of  other  spacecraft  designed  for  near-  earth 


orbital  operations. 


V.  Baseline  Modification  Analysis 


Hie  Modified  Constellation 

Introduction.  The  six-plans,  18-satellite  constellation  analyzed  in  Chapter  IV  was 
recognized  as  having  better  overall  capability  than  a  variety  of  other  candidate  18-satellite 
constellations  evaluated  and  is  now  considered  to  be  the  baseline  OPS  constellation.  One  of  the 
primary  reasons  that  this  constellation  was  chosen  over  an  earlier  three-plane,  nonuniform 
constellation  baseline  proposal  was  that  a  significant  improvement  of  satellite  coverage  was 
obtained  with  the  six-plane  configuration  ( I5-.E9.3.8).  The  99.5*  coverage  of  the  six-plane 
constellation  W8S  found  to  be  significantly  better  than  the  98.4*  coverage  of  the  three-plane 
constellation,  since  the  difference  between  these  two  is  a  measure  of  the  difference  from  total  or 
1 00*  coverage  and  translates  into  size  and  duration  of  degraded  performance  areas  ( 1 5.E9.3.2). 

In  all  satellite  constellations  considered  in  the  selection  process,  circular  orbits  of  equal 
period  were  chosen  as  an  essential  feature;  it  had  been  assumed  that  elliptical  orbits  were 
advantageous  for  coverage  of  limited  areas,  but  that  the  more  uniform  patterns  provided  by 
circular  orbits  were  preferable  for  whole-Earth  coverage  (22:4).  Walker's  expectation  that 
"continuous  whole-Earth  coverage  would  be  provided  most  effectively  by  a  system  in  which  the 
distribution  of  satellites  over  the  Earth's  surface  was  maintained  ss  uniform  as  possible,  subject 
to  the  practical  limitations  imposed  on  a  system  necessarily  involving  multiple  intersecting 
orbits,"  has  been  the  underlying  basis  for  this  assumption.  Although  it  would  appear  to  be  a 
reasonable  assumption,  it  does  not  take  into  account  the  fact  that  system  outages  are  caused  not 
only  by  an  insufficient  number  of  satellites  visible  to  the  user,  but  also  by  the  poor  geometry 
(with  respect  to  the  user)  that  me/  exist  for  those  satellites  in  view. 


Since,  for  the  baseline  constellation,  there  are  always  four  or  more  satellites  in  view  to  the 


earth- based  user ,  all  outages  that  occur  must  be  due  to  the  poor  relative  geometry  of  the  visible 
satellites;  this  was  a  direct  observation  from  the  global  distribution  analysis  of  Chapter  IV.  If  a 
modification  to  this  baseline  constellation  could  be  made  that  would  result  in  more  favorable 
geometry  during  these  outage  periods,  without  causing  an  increase  In  outages  at  other  times,  it 
seems  logical  to  expect  that  coverage  could  be  improved  even  further.  Since  a  change  in  the  shape 
of  the  orbit,  or  eccentricity,  would  result  in  a  change  in  geometry,  it  seems  reasonable  to  expect 
that  elliptical  orbits  might  provide  such  an  opportunity. 

Assumptions.  It  was  assumed  that  the  only  modification  to  be  made  to  the  baseline 
constellation  was  to  the  eccentricity  (shape)  of  the  orbits.  Thus,  all  modified  constellations 
considered  consisted  of  1 8  satellites  deployed  in  six  elliptical  planes,  three  satellites  per  plane, 
each  with  an  orbital  period  of  1 2  hours.  Each  orbital  plane  had  an  inclination  of  55°,  and  was 
separated  from  the  next  by  60  degrees  in  longitude.  The  three  satellites  in  each  orbital  plane  were 
uniformly  distributed  in  time;  this  means  that  each  was  separated  from  the  other  within  a 
particular  orbit  by  four  hours,  or  exactly  one-third  of  the  period.  Relative  phasing  of  the 
satellites  from  one  orbital  plane  to  the  next  remained  at  40°  ;  this  means  that  when  an  ascending 
satellite  in  one  plane  is  crossing  the  equator ,  an  ascending  satellite  in  the  adjacent  plane  to  the  east 
is  40°  above  the  equatorial  plane  in  its  own  orbital  plane.  Since  elliptical  orbits  were  used,  the 
relative  phasing  of  the  satellites  was  approximated  by  positioning  the  ascending  satellite  in  the 
adjacent  plane  to  the  east  ahead  of  the  equator  by  40°  of  eccentric  anomaly-  For  ease  in 
constellation  buildup,  as  well  as  in  the  evaluation  and  comparison  with  the  baseline  constellation, 
it  was  assumed  that  the  eccentricity  chosen  for  one  orbit  would  be  used  for  all  six  orbits. 

Choosing  thfi  Modified  Constellation.  With  the  assumptions  stated  above,  there 
still  remained  the  difficult  task  of  choosing  the  optimum  orbital  eccentricity  and  perigee  locations, 
which,  it  was  hoped,  would  provide  the  improved  coverage  desired.  From  a  practical  standpoint, 
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a  small  eccentricity  would  be  desirable,  as  satellite  stationkeeping  and  antenna  pointing  would  be 
easier  if  the  orbits  were  kept  as  nearly  circular  as  possible.  In  addition,  if  the  orbits  were  made 
too  elliptical ,  outages  would  begin  to  occur  due  to  periods  when  fewer  than  four  satellites  would  be 
visible  to  a  user ,  rather  than  due  solely  to  poor  geometry;  this  would  negate  any  advantage  attained 
by  the  improving  geometry  resulting  from  the  elliptically-shaped  orbits.  Since  this  is  exactly 
what  occurred  at  eccentricities  much  greater  than  0. 1 ,  an  eccentricity  of  .05  was  initially  chosen 
as  a  starting  point  for  the  computer  analysis. 

Since  a  satellite  in  an  elliptical  orbit  travels  at  slowest  speed  near  apogee  and  at  fastest 
speed  near  perigee,  it  spends  the  majority  of  its  time  near  the  apogee  end  of  its  orbit.  By 
positioning  the  apogees  of  each  orbit  over  a  specific  area  of  the  earth  where  increased  coverage  is 
desired,  a  tremendous  improvement  in  performance  was  noted;  unfortunately,  however,  this  was 
at  the  expense  of  a  significant  degradation  of  performance  in  other  arras  of  the  world.  Thus,  when 
the  apogees  of  the  orbits  were  positioned  over  the  middle  northern  latitudes,  the  outages  were 
virtually  eliminated  in  the  northern  hemisphere  while  becoming  more  frequent  and  of  longer 
duration  in  the  southern  hemisphere  (Table  IX).  Asa  result  of  the  preliminary  computer  runs,  it 
was  determined  that  a  near  optimum  location,  on  a  global  basis,  for  the  perigees  of  the  modified 
elliptical  orbits  was  in  the  plane  of  the  equator,  as  outages  occurred  less  frequently  in  these 
particular  cases  than  any  other  location  analyzed. 

After  analyzing  several  data  runs  for  constellations  with  eccentricities  ranging  from  .01  to 
.  t  and  with  all  perigees  located  in  the  equatorial  plane,  an  eccentricity  of  .07  was  selected  as  the 
optimum  value  of  eccentricity  that  would  minimize  the  percentage  of  time  outages  would  occur 
(Table  IX).  It  is  significant  to  note  that  all  eccentricities  between  .0 1  and .  1  provided  superior 
performance  than  that  of  the  baseline  constellation  in  terms  of  reducing  system  outages. 


TABLE  IX 


Global  Distribution  Preliminary  Data  Runs 


Eccentricity 

Arq  of  Perigee  f  Deal 

Hemisphere 

%  PDOP 
>60 

I  PDOP 
>LQ 

.05 

239.7 

Northern 

.0012 

.0010 

.05 

239.7 

Southern 

.0107 

.0087 

.02 

239.7 

Northern 

.0070 

.0048 

.02 

239.7 

Southern 

.0084 

.0068 

.05 

224.5 

Northern 

.0009 

.0003 

.05 

224.5 

Southern 

.0114 

.0106 

.02 

224.5 

Northern 

.0028 

.0023 

.02 

224.5 

Southern 

.0082 

.0069 

.02 

198.4 

Northern 

.0032 

.0025 

.02 

198.4 

Southern 

.0071 

.0062 

.01 

0.0 

Northern 

.0048 

.0039 

.02 

0.0 

Northern 

.0041 

.0037 

.03 

0.0 

Northern 

.0041 

.0034 

.04 

0.0 

Northern 

.0038 

.0029 

.05 

0.0 

Northern 

.0040 

.0026 

.07 

0.0 

Northern 

.0029 

.0017 

.085 

0.0 

Northern 

.0034 

.0015 

.10 

0.0 

Northern 

.0037 

.0016 

Global  Distribution  Run 

Introduction.  Since  the  modified  constellation  selected  was,  like  the  baseline 
constellation,  a  symmetrical  arrangement  of  satellites,  only  the  OOP’s  for  the  northern 
hemisphere  were  calculated.  The  OOP's  for  the  southern  hemisphere  are  identical  for  a  global 
analysis.  The  latitude,  longitude,  and  time  steps  were  chosen  to  be  the  same  as  used  for  the 
baseline  constellation  global  distribution  run,  which  allowed  for  the  direct  comparison  of  the 
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results.  Th8  time  interval  selected  for  the  run,  as  in  the  case  of  the  baseline  constellation,  was 
four  hours,  as  after  this  length  of  time  the  pattern  of  outages  was  repetitive  due  to  the  symmetric 
arrangement  of  satellites. 


TABLE  X 

Modified  Constellation  Orbital  Elements 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


ORBITAL  ELEMENTS 

ECC  AR6P  RASC  INC  ANOM  PER 
(Deg)  (Deg)  (Deg)  (Deg)  (Hrs) 


0.07 

0.00 

30.00 

0.07 

0.00 

30.00 

0.07 

0.00 

30.00 

0.07 

0.00 

90.00 

0.07 

0.00 

90.00 

0.07 

0.00 

90.00 

0.07 

0.00 

150.00 

0.07 

0.00 

150.00 

0.07 

0.00 

150.00 

0.07 

0.00 

210.00 

0.07 

0.00 

210.00 

0.07 

0.00 

210.00 

0.07 

0.00 

270.00 

0.07 

0.00 

270.00 

0.07 

0.00 

270.00 

0.07 

0.00 

330.00 

0.07 

0.00 

330.00 

0.07 

0.00 

330.00 

55.00 

0.00 

12.00 

55.00 

126.60 

12.00 

55.00 

233.40 

12.00 

55.00 

42.80 

12.00 

55.00 

160.30 

12.00 

55.00 

269.50 

12.00 

55.00 

83.80 

12.00 

55.00 

193.90 

12.00 

55.00 

310.00 

12.00 

55.00 

123.30 

12.00 

55.00 

230.00 

12.00 

55.00 

355.70 

12.00 

55.00 

161.40 

12.00 

55.00 

270.60 

12.00 

55.00 

44.10 

12.00 

55.00 

198.60 

12.00 

55.00 

315.90 

12.00 

55.00 

89.40 

12.00 

PARAMETERS  USED  IN  GLOBAL  DISTRIBUTION  CALCULATIONS: 

MASKING  ANGLE  «  5.00  DEGREES 
LATITUDE  STEP  -  10.00  DEGREES 
LONGITUOE  STEP  -  20.00  DEGREES 
TOTAL  TIME  (MIN)  -  240 
TIME  INCREMENT  (MIN)  -  5 
HEMISPHERE  EVALUATED  -  NORTHERN 
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Results.  The  results  of  the  global  analysis  computer  run  for  the  modified  constellation 
are  shown  in  Table  XI.  It  is  significant  to  note  that,  for  the  modified  constellation,  outages  occur 
only  0.29 *  of  the  time,  or  stated  in  another  way,  PDOP  values  of  less  than  6.0  are  available 
99.71  *  of  the  time;  this  represents  a  47*  reduction  of  outages  over  the  baseline  constellation 
(Figure  1 4),  which  h8S  outages  occuring  0.55*  of  the  time.  If  PDOP  values  greater  than  7.0  are 
used  as  a  means  of  comparison,  outages  ( PDOP  values  greater  that  7.0)  are  reduced  by  over  56*. 


TABLE  XI 

Modified  Constellation  Global  Distribution  -  DOP  Values 


DILUTION  OF  PRECISION  PARAMETERS  -  ACCUMULATIVE  GLOBAL  DISTRIBUTION 
( Percentage  of  T ime  That  DOP  Value  Shown  is  Exceeded) 


■!  »!•  :■§ 

■a  iTTiH 

mw-m 

u^m 

MW* 

.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

.0000 

1.0000 

.7297 

1.0000 

.5601 

.0439 

.1009 

.9849 

.9991 

.0030 

.0018 

.0143 

.2257 

.0334 

.0009 

.0018 

.0650 

.1111 

.0063 

.0004 

.0004 

.0013 

.0138 

.0431 

.0022 

.0004 

.0003 

.0011 

.0029 

0114 

.0012 

.0003 

.0017 

.0028 

The  highlighted  number  represents  the  percentage  of  time  that  an  outage  occurs  due 
to  poor  satellite  geometry  of  the  baseline  constellation.  A  PDOP  value  greater  than  6.0  is 
considered  to  constitute  such  an  outage. 


This  reduction  in  system  outages  that  occurs  with  the  modified  system  is  not  without  some 
cost.  The  cost  of  this  increased  coverage,  however,  is  limited  to  a  very  slight  reduction  in  the 
percentage  of  times  that  the  best  POOP  values  are  obtainable,  as  depicted  in  Figure  15.  An 
examination  of  the  data  in  Table  XII  and  Table  Xltt  provides  the  reason  for  this  slight  degradation  of 
the  best  POOP  values  obtainable.  Although  the  minimum  and  maximum  number  of  satellites 
visible  to  the  global  user  remained  basically  the  same  as  for  the  baseline  constellation,  the 
distribution  of  their  occurrences  did  not.  Since  the  analysis  showed  that  there  were  always  four 
or  more  satellites  visible  to  the  earth-based  user  for  both  the  baseline  constellation  and  the 
modified  system,  the  outages  that  occurred  were  attributed  to  the  poor  relative  geometry  of  the 
four  satellites  selected.  In  almost  every  case  where  this  occurred  for  the  baseline  constellation, 
there  were  only  the  minimum  number  of  four  satellites  from  which  to  choose  from.  A  close 
examination  of  Table  Yl  confirms  this,  as  the  data  clearly  shows  that  five  or  more  satellites  were 
in  view  1 00X  of  the  time  at  all  northern  latitudes  sampled  except  at  60°N  and  40°N,  which  were 
also  the  latitudes  at  which  the  outages  occurred. 

The  modified  constellation,  however,  provides  a  higher  probability  of  seeing  five  or  more 
satellites  at  these  same  latitudes,  which  means  that  in  addition  to  the  change  in  relative  geometry 
of  satellites  provided  by  the  elliptical  orbits,  additional  satellites  from  which  to  choose  are  made 
available  over  the  specific  regions  of  the  earth  where  the  outages  had  occurred.  Consequently,  the 
addition  of  a  fifth  satellite  over  these  areas  provided  an  opportunity  for  a  better  relative  geometry 
of  the  four  satellites  selected  and  eliminated  nearly  half  the  outages  that  had  previously  occurred 
in  these  areas.  Since  this  fifth  satellite  was  provided  at  the  expense  of  a  reduction  of  the  number 
of  satellites  available  in  other  areas  which  previously  had  a  higher  number  of  satellites  from 
which  to  select,  the  best  POOP  values  obtainable  in  these  other  areas  were  slightly  degraded. 
Thus,  although  on  a  global  basis  the  modified  constellation  provided  a  slightly  lower  probability  in 


51 


percent  (98.92*  to  99.17*)  that  five  or  more  satellites  would  be  seen,  it  effectively 
redistributed  the  "wealth"  by  targeting  those  areas  of  the  world  with  weakest  coverage. 


TABLE  XII 

Maximum  and  Minimum  Numbers  Seen  at  Each  Latitude  and  Longitude- Modified  Constellation 


MAXIMUM 
LONGITUDE  (DEG) 

1  1  1  1  122222333 

24680246802468024 
LAT  000000000000000000 

90  888888888888888888 

80  788878887887788887 

70  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8  8 

60  888888888888888888 

50  777777776777677777 

40  7  6  6  8  7  7  6  8  8  7  7  6  8  7  6  6  8  8 

30  7  7  6  7  7  7  6  8  8  7  7  6  7  7  7  6  7  7 

20  687768877777768877 

10  787888877787778877 

0  787787778787878877 


Probability  of  Seeing  N  or  More  Satellites-Modified  Constellation 


NUMBER  OF  SATELLITES 


LAT 

0 

1 

2 

3 

4 

5 

6 

7 

8 

PROBABILITY  (IN  PERCENT) 

OF  SEEING  N  OR  MORE  SATELLITES 

90 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

40.82 

24/ 

80 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 
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ON  A  6L06AL  BASIS  THE  PROBABILITY  (IN  PERCENT)  THAT  N  OR  MORE  SATELLITES  WILL  BE  VISIBLE 

NUMBER  OF  SATELLITES 


System  Degradation  Qj]£  Satellite  Losses 

Introduction  The  effect  of  satellite  losses  from  the  modified  baseline  constellation  was 
analyzed  for  comparison  with  the  results  obtained  from  the  analysis  of  the  baseline  constellation 
discussed  in  Chapter  IV.  For  this  limited  analysis,  the  effect  upon  geometric  performance  from 
losses  of  one  to  three  satellites  was  evaluated  for  the  same  cases  that  provided  both  the  best  and 
worst  performance  for  the  baseline  constellation.  Whether  or  not  these  cases  are  the  equivalent 
best  and  worst  cases  for  the  modified  constellation  as  well  remains  to  be  determined,  but  since  this 
part  of  the  analysis  was  conducted  solely  for  the  purpose  of  comparing  the  two  constellations,  these 
cases  should  provide  adequate  representation  of  the  broad  range  of  effects. 

For  each  computer  run,  a  ten  minute  time  increment  was  selected  and  the  constellation  was 
evaluated  over  a  six  hour  period.  The  latitude  and  longitude  step  sizes  were  20°  and  both 
hemispheres  were  evaluated.  These  parameters  were  identical  to  those  used  in  the  analysis  of 
satellite  losses  from  the  baseline  constellation ,  which  adds  validity  to  the  comparison. 

Results.  The  percentage  of  time  that  system  outages  occurred  due  to  satellite  losses  from 
the  modified  constellation  are  listed  in  Table  XIV.  The  percentages  in  brackets  represent  the 
corresponding  values  obtained  for  the  same  losses  analyzed  for  the  baseline  constellation  (Table 
VI I ).  In  practically  every  case,  the  modified  baseline  constellation  provided  superior  performance 
over  the  baseline  constellation.  As  seen  in  the  baseline  analysis,  the  particular  satellites  removed 
from  the  constellation  made  a  large  difference  in  the  amount  of  degradation  occurring;  outages 
occurred  in  the  three  satellite  loss  case  from  8. 11 51  to  15.3*  of  the  time  and  in  the  two  satellite 


case  from  4.555*  to  7.47*  of  the  time,  representing  the  best  and  worst  case  scenarios. 


TABLE  XIV 


System  Performance  Degradation  Due  to  Satellite  Losses-  Modified  Constellation 


Satellites 

%  PDOP 

2  PDOP 

Deleted 

>  6.Q 

LU1 

1 

.0228  [.0232] 

.0161  [.0201] 

1.5* 

.0455  [.0478] 

.0329  [.0407] 

1,10** 

.0774  [.0800] 

.0661  [0725] 

4,5,6* 

.0811  [.0780] 

.0574  [.0679] 

1,8,15** 

.1530  [.1584] 

.1349  [1454] 

*  Best  Case 
**  Worst  Case 


Satellite  Usee 

Introduction.  In  analyzing  the  geometric  performance  of  the  modified  baseline 
constellation  for  the  satellite  user,  three  cases  were  considered:  an  intermediate  altitude  earth 
orbit,  a  high  altitude  earth  orbit,  and  a  high  altitude  ballistic  missile  trajectory.  The  low  altitude 
earth  orbit  and  ICBM  trajectory  were  not  considered  for  this  analysis  as  it  was  assumed  that  the 
modified  constellation  would  provide  outstanding  geometric  performance  for  these  two  users,  as 
had  been  the  case  with  the  baseline  constellation.  In  any  event,  the  three  cases  selected  were 
representative  of  the  wide  range  of  potential  satellite  users  and  were  considered  adequate  for  the 
purposes  of  this  limited  comparison. 

The  parameters  selected  for  the  computer  analyses  of  each  case  were  identical  with  those 
used  in  the  analyses  of  these  cases  for  the  baseline  constellation,  and  the  cases  analyzed  were 
identical  to  those  described  in  Chapter  IV.  It  was  assumed  that  the  OPS  antennas  remained  pointed 


at  the  center  of  the  earth  throughout  the  orbit  of  the  satellites,  and  that  a  sufficiently  strong  signal 
was  available  to  the  user  satellite  at  all  distances. 


Intermediate  Altitude  Earth  Orbits.  The  orbital  parameters  of  the  highly 
elliptical  satellite  user  orbit  considered  were  described  in  Chapter  IV  (Case  2).  It  was  analyzed 
for  an  antenna  half-angle  beamwidth  of  21.4  degrees  and  45  degrees.  With  the  designed  antenna 
half-angle  beamwidth  of  21.4  degrees,  four  or  more  satellites  were  visible  18.6*  of  the  time 
compared  with  1 7.9*  for  the  baseline  constellation.  Three  or  more  satellites  were  visible  26.2* 
of  the  time  versus  24. 1  *  for  the  baseline,  and  two  or  more  satellites  were  visible  53.8*  of  the 
time  compared  with  5 1  *  for  the  baseline.  No  satellites  were  visible  to  the  satellite  user  27.6* 
of  the  time. 

When  the  same  user  orbit  was  analyzed  using  an  antenna  half-angle  beamwidth  of  45 
degrees,  as  in  the  case  of  the  baseline  analysis,  significantly  better  results  were  obtained.  Four  or 
more  satellites  were  available  to  the  user  at  all  times  during  its  orbit,  with  POOP  values  ranging 
from  1.53  to  123.  PDOP  values  less  than  6.0  were  achieved  32.4*  of  the  time  compared  with 
35.2*  for  the  baseline  constellation.  The  results  are  summarized  in  Table  XV. 

High  Altitude  Ballistic  Missile  Trajectory.  The  same  highly  elliptical  missile 
trajectory  (Case  4)  described  in  Chapter  IV  was  used  in  this  analysis  as  had  been  used  for  the 
baseline  analysis.  Half-angle  beamwidths  of  21.4  and  45  degrees  were  selected  for  the 
comparison  evaluation.  With  the  designed  antenna  half-angle  beamwidth,  four  or  more  satellites 
were  visible  to  the  user  7.8*  of  the  time,  identical  to  the  results  of  the  baseline  analysis.  Three 
or  more  satellites  were  available  1 4.0*  of  the  time  compared  to  the  baseline  1 2.4* ,  two  or  more 
were  visible  30.2*  of  the  time  versus  29.5*  for  the  baseline,  and  one  or  more  satellites  were 
available  61.2*  of  the  time  compared  to  62.8*  for  the  baseline  constellation.  There  were  no 
satellites  visible  38.8*  of  the  time.  When  the  half-angle  beamwidth  was  increased  to  45  degrees, 


four  or  more  satellites  were  always  visible,  as  was  found  to  be  the  case  for  the  baseline 
constellation.  PDOP  values  less  than  6.0  were  achieved  in  this  case  22.5$  of  the  time,  which  was 
a  slight  improvement  over  the  20.9$  achieved  by  the  baseline  constellation. 


Table  XV 

Modified  Constellation  Evaluation  -  Satellite  User 


$  of  Time  n  or  More  Satellites 


User 

Antenna 

Visible 

Satellite 

BWIDTH 

1 

2 

3 

4 

PDOP  Range 

Case  2 

21.4 

72.4 

53.8 

26.2 

18.6 

1.536- 

45.0 

100.0 

100.0 

100.0 

100.0 

1.539-  123 

Case  4 

21.4 

61.2 

30.2 

14.0 

7.8 

1.506- 

45.0 

100.0 

100.0 

100.0 

100.0 

1.598-345 

Case  5 

21.4 

53.8 

25.5 

2.7 

0.0 

45.0 

100.0 

100.0 

100.0 

100.0 

14.4-  166.2 

90.0 

100.0 

100.0 

100.0 

100.0 

5.7  -  8.8 

High  Altitude  Earth  Orbits.  A  circular  geosynchronous  orbit  W8S selected  (Case  5) 
for  evaluation  at  antenna  half-angle  beamwidths  of  21.4, 45,  and  90  degrees.  With  the  designed 
antenna  beamwidth,  there  were  never  four  satellites  visible.  Three  or  more  satellites  were 
visible  2.7$  of  the  time  for  the  modified  constellation  compared  with  1.4$  of  the  time  fur  the 
baseline.  Two  or  more  satellites  were  visible  25.5$  of  the  time  and  one  or  more  visible  only 
53.8$  of  the  time,  compared  with  22. 1  $  and  55.2$  respectively  for  the  baseline. 


When  the  antenna  half-angle  beamwidth  was  increased  to  45  degrees  and  higher,  four  or 
more  satellites  were  always  available  to  the  user.  At  a  half-angle  beamwidth  of  90  degrees,  PDOP 
values  less  than  7.0  were  achieved  75.8*  of  the  time  compared  with  62. 1  *  of  the  time  for  the 


baseline  constellation. 

Summary.  The  geometric  performance  of  the  modified  constellation  for  the  satellite  user 
is  probably  slightly  favorable  to  that  of  the  baseline  constellation.  For  all  three  cases  analyzed, 
the  modified  constellation  generally  provides  more  satellites  available  to  the  user  when  less  than 
four  satellites  are  available.  When  more  than  four  satellites  are  visible  to  the  user,  PDOP  values 
less  than  6.0  or  7.0  are  achieved  a  higher  percentage  of  time  than  for  the  unmodified  baseline 
constellation,  resulting  in  fewer  outages.  Although  only  a  limited  number  of  representative  orbits 
were  considered  in  the  analysis,  one  can  only  conclude  from  the  results  that  the  modified 
constellation  performs  at  least  on  an  equal  basis  with  that  of  the  baseline  constellation  for  the 
typical  satellite  user,  and  probably  better.  Although  GPS  performance  for  the  space-based  user 
was  not  considered  as  a  factor  in  the  selection  of  the  baseline  constellation,  it  certainly  merits 
consideration.  Both  the  baseline  constellation  and  the  modified  constellation  show  great  potential 
for  use  in  near-earth  space  navigation. 


VI.  Conclusions/Recommendations 

The  Baseline  Satellite  Constellation 

The  1 8-satellite  configuration  now  considered  the  GPS  baseline  constellation  will  provide 
nearly  continuous  coverage  on  a  worldwide  basis  to  its  earth- based  users,  providing  them  with 
positioning  capability  of  high  accuracy  on  the  order  of  tens  of  meters.  Those  system  outages  that  do 
occur  with  the  proposed  system  will  be  due  solely  to  the  poor  geometry  of  the  selected  satellites, 
and  these  outages  will  only  occur  (on  an  average)  approximately  one-half  percent  of  the  time.  The 
duration  of  these  outages  will  be  brief  ( 5-30  minutes)  and  the  locations  predictable. 

Satellite  losses  sustained  by  the  18-satellite  constellation  will  significantly  degrade  the 
coverage  available,  and  the  amount  of  degradation  caused  by  such  losses  will  be  largely  dependent 
on  the  relative  positions  of  those  satellites  lost.  Should  replacement  spares  be  unavailable, 
rephasing  the  remaining  satellites  could  substantially  improve  the  coverage  until  the  constellation 
is  restored  to  its  original  number.  The  planned  placement  of  three  active,  in-orbit  spares  should 
greatly  increase  the  system’s  reliability  and  provide  some  flexibility  in  rephasing  satellites  to 
compensate  for  unplanned  losses. 

For  the  satellite  or  space- based  user,  the  feasibility  of  using  GPS  for  positioning  in  the 
conventional  way  is  primarily  a  function  of  both  the  user's  altitude  and  the  antenna  beamwidth  of 
the  GPS  satellites.  With  the  currently  designed  antenna  half-angle  beamwidth  of  2 1 A  degrees,  the 
low-  altitude  satellite  user  will  usually  have  positioning  accuracy  capability  superior  to  that  of 
the  typical  global  user  due  to  the  incre&>jd  number  of  satellites  from  which  to  select  and  the 
resulting  better  geometry.  Due  to  the  limiting  constraints  of  the  antenna  beamwidth ,  however ,  the 
space-based  user's  positioning  accuracy  will  rapidly  deteriorate  with  increasing  altitude,  as 
fewer  and  fewer  satellites  remain  available  for  selection.  Outages  become  more  and  more  frequent 
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until  altitudes  are  reached  where  conventional  positioning  using  four  satellites  simultaneously 
becomes  impossible.  Modifying  the  antenna  beamwidth  design  of  the  OPS  satellites  appears  very 
promising  for  expanding  the  usefulness  of  OPS  to  the  high  altitude  satellite  user ,  but  even  without 
such  a  modification,  greatly  improved  positioning  capability  for  these  users  can  be  obtained  by 
talcing  sequential  measurements  from  OPS  satellites  as  the/  become  available.  Equipping  the 
receiver  with  a  precise  clock  and  incorporating  these  sequential  measurements  with  the  memory 
inherent  in  the  known  orbital  dynamics  of  the  satellite  will  make  navigation  possible  for  these 
users  as  well. 


The  Baseline  Modification 

In  an  attempt  to  improve  the  geometric  performance  of  the  baseline  constellation  by 
reducing  the  number  and  duration  of  system  outages,  a  modification  was  made  to  the  baseline 
constellation  which  changed  the  shape  of  the  circular  GPS  satellite  orbits  to  slightly  elliptical 
ones,  each  with  an  eccentricity  of  .07.  The  perigee  of  each  satellite  orbit  was  positioned  in  the 
plane  of  the  equator;  all  other  orbital  parameters  remained  unchanged.  This  modification  reduced 
the  percentage  of  system  outages  occurring  on  a  global  basis  by  nearly  502,  which  was  a 
significant  improvement  in  coverage  over  the  baseline  constellation.  The  "cost"  for  this 
improvement  was  a  slight  reduction  in  the  best  accuracies  attainable  for  sane  areas  of  the  world. 

Although  elliptical  orbits  are  certainly  advantageous  for  coverage  of  limited  areas  of  the 
world,  it  has  generally  been  assumed  that  circular  orbits  are  preferable  for  whole-earth 
coverage.  This  assumption,  however ,  does  not  appear  to  be  valid  at  leest  in  this  particular  case,  os 
the  outages  that  occurred  with  the  baseline  constellation  were  due  solely  to  the  poor  geometry  of 
the  available  satellites,  and  not  to  the  lack  of  sufficient  satellites  for  measurements.  Changing  the 
eccentricity  of  the  satellite  orbits  resulted  in  a  more  favorable  satellite  geometry  relative  to  the 


users  on  a  global  average;  consequently,  fewer  outages  occurred. 

The  modified  constellation  also  compared  favorably  to  the  baseline  constellation  with  respect 
to  the  space-based  user,  providing  positional  accuracy  generally  better  than  that  provided  to  the 
user  by  the  unmodified  constellation.  Since  only  a  limited  number  of  user  orbital  trajectories 
were  examined,  it  can  only  be  concluded  that  the  modified  constellation  provides  the  space-based 
user  with  positioning  accuracy  comparable  to  that  of  the  baseline,  and  perhaps  better . 

The  effect  of  satellite  losses  from  the  modified  constellation  were  also  examined  for  the  cases 
which  provided  both  the  best  and  worst  navigational  performance  for  the  baseline  constellation.  In 
almost  every  case,  the  modified  constellation  provided  better  overall  coverage  than  the  baseline 
constellation. 

Recommendations 

Although  the  location  and  duration  of  system  outages  is  considered  a  significant  factor  in 
selecting  the  optimum  18-satellite  OPS  constellation,  other  factors  may  make  the  proposed 
modification  to  the  baseline  constellation  impractical  to  implement.  It  has  been  assumed  that  the 
necessary  satellite  station-keeping  is  available  to  maintain  the  orbital  pattern  for  the  expected 
lifetime  of  the  satellite;  precession  of  the  line  of  apsides  (perigee  movement)  due  to  earth 
oblateness,  for  example,  is  approximately  .022  degrees  per  day,  and  would  require  periodic 
corrections.  Further  analysis  is  recommended  to  determine  the  feasibility  of  this  modification  to 
the  proposed  baseline  constellation.  Since  the  primary  intent  of  this  analysis  wss  to  determine 
whether  changing  the  eccentricity  of  the  OPS  satellite  orbits  could  improve  geometric 
performance  and  reduce  the  percentage  of  system  outages,  and  not  to  find  the  optimal 
constellation,  other  elliptical  variations  of  the  baseline  should  also  be  analyzed  for  comparison. 

The  feasibility  of  modifying  OPS  antenna  design  to  better  accommodate  the  space- based  user 
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should  also  be  closely  examined.  Although  GPS  was  not  originally  designed  for  use  in  space 
navigation,  the  potential  applications  into  this  area  appear  very  promising.  Increasing  the 
antenna  beamwidth,  if  practical,  would  greatly  improve  the  availability  of  the  GPS  satellites  to 
users  operating  in  space,  and  would  make  available  to  the  user  extremely  accurate  estimates  of 
position 
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Appendix  A 

Glossary  of  Technical  Terms 


Antenna  Half- angle  Beamwldth  -  the  angle  between  the  antenna  centerline  and  the  edge  of  the 
main  lobe  skirt  of  the  OPS  antenna  pattern 

Apogee  -  the  point  in  a  satellite's  orbit  farthest  from  the  central  attracting  body  (earth) 

Argument  gf  Perioee  -  the  angle,  in  the  plane  of  the  satellite's  orbit,  between  the  ascending  node 
and  the  periapsis  point  (perigee),  measured  in  the  direction  of  satellite  motion 

Ascending  Node  -  the  vector  pointing  from  the  center  of  the  earth  to  the  satellite  as  it  passes 
through  the  fundamental  (or  equatorial)  plane  in  a  northerly  direction 

Clock  Bias  -  the  error,  or  fixed  bias,  in  the  user's  imprecise  clock  relative  to  OPS  system  time 
provided  by  the  OPS  satellites  which  contributes  an  error  in  range  measurement 

Dilution  fll  Precision  (POP)  -  a  parameter  used  for  measuring  system  position  accuracy  of 
the  OPS  satellites  as  a  function  of  the  variability  of  satellite  geometry 

Direction  Cosines  -  the  numbers,  representing  the  cosines  of  three  angles,  which  completely 
define  the  direction  of  a  given  vector  relative  to  a  given  coordinate  system 

Eccentric  Anomaly  -  the  angle  between  the  semi-major  axis  of  an  ellipse  and  a  point  on  a  circle 
circumscribed  about  the  ellipse  represented  by  the  interesection  of  that  circle  with  a  line, 
perpendicular  to  the  semi-major  axis,  through  the  point  of  interest  on  the  ellipse 

Eccentricity  -  a  constant  defining  the  shape  of  an  elliptical  or  other  conic  orbit 

Fundamental  Plane  -  the  plane  of  the  equator 

800P  -  the  geometric  dilution  of  precision,  a  parameter  which  reflects  the  dilution  of  precision  of 
position  accuracy  in  three  dimensions  plus  time 

Oaoayncfironous  Orbit  -  a  circular  orbit  with  a  period  of  24  hours,  situated  in  the  plane  of  the 
equator ,  such  that  the  satellite  always  remains  over  the  same  point  of  the  earth 

Hohman  Transfer  -  an  elliptical  transfer  orbit  between  two  circular  orbits  requiring  minimum 
fuel  for  the  maneuver 

HDOP  -  the  horizontal  dilution  of  precision,  a  parameter  which  reflects  the  dilution  of  precision  of 
position  accuracy  in  the  two  horizontal  directions 


Inclination  -  the  angle  which  describes  the  orientation  of  the  orbital  plane  with  the  plane  of  the 
equator 

Ionospheric  Delay  -  the  time  delay  of  RF  signals  passing  through  the  ionosphere  due  to  a 
reduction  of  speed  and  bending  of  the  ray  from  refraction,  Inversely  proportional  to  the  square  of 
the  frequency 

Julian  Tima  -  the  time  from  which  the  Julian  Calendar  dates  ( 46  B.C. ) 

Longitude  o£  Ascending  Node  -  the  angle,  in  the  plane  of  the  equator,  between  the  vector 
pointing  to  the  vernal  equinox  and  the  ascending  node,  measured  counterclockwise  when  viewed 
from  the  north  side  of  the  fundamental  plane 

Hashing  Angle  -  the  minimum  elevation  angle,  due  to  terrain  obstructions  and  the  earth's 
atmosphere,  that  a  GPS  satellite  can  have  relative  to  the  user  and  still  be  usable 

MDOP  -  a  parameter  which  reflects  the  dilution  of  precision  of  position  accuracy  in  the  larger 
component  of  the  horizontal  position  error 

Orbital  Elements  -  the  six,  independent  quantities  which  are  sufficient  to  completely  describe 
the  size,  shape,  and  orientation  of  on  orbit  and  the  position  of  the  satellite  along  the  orbit  at  a 
particular  time 

Outage  -  a  situation  occurring  for  a  user  at  a  specific  location  and  time  period  when  positioning 
capability  is  unavailable  due  to  either  unfavorable  satellite  geometry  or  an  insufficient  number  of 
visible  GPS  satellites  ( less  than  four) 

PDOP  -  the  position  dilution  of  precision  parameter  which  reflects  the  dilution  of  positioning 
accuracy  in  three  dimensions 

Perigee  -  the  point  in  a  satellite's  orbit  closest  to  the  earth 

Period  -  the  time  required  for  a  satellite  to  traverse  its  entire  orbit  one  time,  dependent  only  on 
the  size  of  the  semi-major  axis 

Pseudoranoe  -  the  sum  of  the  actual  range  displacement  from  the  user  to  GPS  satellite  plus  the 
offset  due  to  the  user  time  error  (clock  bias) 

Relotlve  Phasing  -  the  number  of  degrees  that  an  ascending  satellite  in  one  orbital  plane  is  above 
the  equatorial  plane  relative  to  an  ascending  satellite  which  is  crossing  the  equatorial  plane  in  an 
adjacent  orbital  plane  to  the  west 

Satellite  Constellation  -  a  particular  arrangement  or  configuration  of  a  specific  number  of 
satellites 

Satellite  Geometry  -  the  geometrical  relationship  of  the  four  selected  GPS  satellites  relative  to 
the  user  and  each  other 


Satellite  Visibility  -  a  measurement  of  the  number  of  satellites  above  a  minimum  elevation 
angle  ( masking  angle)  visible  to  a  user  at  a  specific  time  and  location 

Seml-malor  Axis  -  a  constant  defining  the  size  of  a  conic  orbit 

Satellite  User  -  a  OPS  user  situated  in  space  (space- based  user) 

Sidereal  Time  -  the  time  based  on  a  sidereal  day,  which  is  defined  as  the  time  required  for  the 
earth  to  rotate  once  on  its  axis  relative  to  the  stars  ( approximately  23h56m04s  of  solar  time) 

TPOP  -  the  dilution  of  precision  in  time  which  is  an  estimate  of  the  range  equivalent  of  the  user 
clock  bias 

True  Anomaly  -  the  angle,  in  the  plane  of  the  satellite's  orbit,  between  periapsis  and  the  position 
of  the  satellite  at  a  particular  time,  or  epoch 

Universal  Time  -  the  local  mean  solar  time  on  the  Greenwich  meridian,  also  called  Greenwich 
Mean  Time  or  Zulu  time 


YPOP  -  the  parameter  reflecting  the  dilution  of  precision  in  the  vertical  dimension  (altitude 


Appendix  B 

The  Computer  Program 


Introduction 

The  computer  program  used  for  the  analysis  of  all  satellite  constellations  evaluated  in  this 
stud/  is  a  modification  of  a  Fortran  program ,  developed  by  the  Rand  Corporation,  on  the  geometric 
performance  of  pseudoranging  satellite  systems  ( 1 7).  Since  a  complete,  detailed  description  of 
the  original  program  and  its  operation  can  be  found  in  Reference  1 7,  the  information  contained  in 
this  appendix  will  be  limited  to  an  explanation  of  the  variables  aid  subroutines  used  in  the 
program  and  the  modifications  made.  In  addition  to  providing  the  complete  program  listing  in 
Appendix  C,  a  sample  of  the  output  for  the  two  types  of  cases  (global  distribution  and  satellite 
user  calculations)  analyzed  in  this  study  is  provided  in  Appendix  D. 

There  are  three  types  of  calculations  performed  by  the  program:  in  Case  I ,  the  user  is  on  a 
satellite;  in  Case  II,  the  user  is  positioned  at  a  specific  latitude  and  longitude  on  the  earths 
surface;  and,  in  Case  III ,  a  group  of  users  are  located  at  a  set  of  latitudes  and  longitudes  forming  a 
net  over  the  whole  surface  or  an  entire  hemisphere  of  the  earth  ( 1 7: 1 5).  Since  only  Case  I  and 
Case  III  were  used  in  this  evaluation,  no  sample  output  is  provided  for  Case  II;  that  portion  of  the 
program  was  not  modified  and  will  not  be  addressed  in  this  discussion. 

Program  Modifications 

The  original  program  was  written  in  Fortran  IV  and  implemented  on  an  IBM  370/158 
computer.  This  version  of  the  program  has  been  converted  to  Fortran  V  for  compatibility  with  a 
Fortran  V  compiler.  The  size  dimension  of  matrix  KXX  in  the  main  program  and  in  subroutines 
TAAT  and  TALL  was  reduced  to  accomodate  a  maximum  of  24  navsats,  rather  than  the  36  navsats 


provided  for  in  the  original  program ;  this  was  done  to  cut  down  on  the  memory  core  storage  used  in 
the  program's  operation.  A  local  library  matrix  inversion  routine  (LINV2F),  was  substituted  for 
the  matrix  inversion  routine  used  in  subroutine  COVNAV  of  the  original  program ,  and  for  greater 
ease  in  entering  data,  formatting  was  changed  to  allow  free  field  data  input. 

Although  the  original  program  contained  statements  that  would  allow  the  program  to 
continue  operation  when  less  than  four  satellites  were  available  to  the  earth- based  user,  it  did  not 
consider  these  occurrences  in  the  calculations  of  OOP  values  and  thus  provided  erroneous  data  for 
the  global  distribution  calculations  when  this  occurred.  The  program  was  modified  to  incorporate 
these  situations  into  the  probability  calculations  of  OOP  values  and  thus  provide  accurate  data  For 
the  satellite  user ,  it  was  discovered  that  the  original  program's  output  was  not  affected  by  a  change 
in  antenna  beam  width;  this  problem  was  corrected  in  the  current  version.  A  modification  was  also 
made  to  the  satellite  user  portion  of  the  program  to  provide  as  output  the  number  of  navsats 
available  to  the  user  when  less  than  four  satellites  were  visible;  the  original  program  did  not 
provide  this  information  as  output 

Since  average  OOP  values  were  not  calculated  by  the  program  for  the  global  distribution 
runs  (the  primary  interest  was  in  determining  when  the  DOP  values  would  exceed  certain  limits 
and  cause  system  outages  to  occur),  additional  statements  were  added  to  the  program  setting  default 
DOP  values  equal  to  1 0OO  when  the  covariance  matrix  became  singular  due  to  poor  geometry.  This 
allowed  the  program  to  continue  operation  when  such  situations  were  encountered,  while  the 
original  program  would  have  terminated.  Since  ,  in  reality,  the  DOP  values  approach  infinity  as 
an  outage  occurs,  and  DOP  values  much  less  than  this  are  considered  as  constituting  a  system 
outage  ( POOP  values  less  than  6,  for  example)  ,  the  exact  values  are  unimportant  at  these  points 
as  long  as  they  are  set  substantially  higher  than  those  values  constituting  such  an  outage.  The 
EGAD  program  currently  in  use  by  the  Aerospace  Corporation  employs  the  same  technique. 


One  additional  subroutine  (SATDEL)  was  added  to  the  program  to  facilitate  the  analysis  of  the 
effect  of  satellite  losses  on  a  given  constellation.  Subroutine  SATDEL  allows  for  the  program  user 
to  selectively  delete  from  one  to  ten  satellites  from  the  original  constellation  for  use  in  comparing 
the  effect  upon  geometric  performance  caused  by  a  particular  combination  of  satellites.  The 
number  of  satellites  to  be  deleted  and  the  particular  identification  numbers  of  these  satellites  are 
entered  as  part  of  the  data  input  at  the  beginning  of  the  program.  This  subroutine  was  necessary  to 
identify  which  combination  of  satellites  provided  the  best  and  worst  case  navigational  accuracy. 

Explanation  of  Variables  In  Main  Program 

The  variables  used  in  the  main  program  ami  an  explanation  of  their  use  is  provided  in  Table 
XVI  and  Table  XVI  I.  Table  XYI  lists  and  explains  the  variables  used  for  that  portion  of  the  program 
involving  the  space- based  or  satellite  user ,  while  Table  XVI I  explains  those  variables  used  for  the 
global  distribution  calculations  involving  the  user  on  the  ground. 


TABLE  XVI 


Explanation  of  Variables  in  Main  Program-  Satellite  User  ( 1 7: 19-20) 

FORTRAN 

EQUATIONS 

EXPLANATION 

RF(IV) 

RMX(N.IV) 

R(IY) 

£ 

Vector  from  the  center  of  the  earth  to  a  navsat  ( where 
N  is  the  "identification  number"  of  the  navsat  and  "IV" 
is  the  index  of  the  components  of  the  vector) 

UPV(IV) 

RMX(IP.IV) 

? 

Vector  from  the  center  of  the  earth  to  a  user  satellite 
(where  “IP"  is  the  identification  number"  of  the  user 
satellite  and  "IV"  is  as  above) 

UTS 

0 

Vector  from  a  user  satellite  to  a  navsat 

AR 

R 

Length  of  vector  £ 

AP 

p 

Length  of  vector  P 

A 

Cos-1(r/R),  where  r  =  radius  of  earth 

B 

Cos-,(r/P) 

PHI 

♦ 

Cos'Ur/R)  +  cos^Ur/P) 

THETN 

*n 

Cos-,(£»£/RP) 

AU 

U 

Length  of  vector  0 

BETAN 

Pn 

Cos_,(-0«£/UP) 

AIWID 

AIN 

<x 

Navsat  antenna  beemwidth  half- angle  ( input)  relative 
to  a  vector  from  the  navsat  to  the  center  of  the  earth 

TABLE  XVI  (Continued) 


FORTRAN  EQUATIONS  EXPLANATION 


BWIDTH 

^  “  ^®n~  fV 

USUV(IINT(N),IV) 

Bj,  i=1 ,2,3,4 

Unit  vectors  from  a  user  toward  4  navsats  (where 

U(IV) 

IINT(N)  contains  the  "identification  number"  of  the 
satellites  which  have  not  been  eliminated  and  "IV 
as  before) 

THETT 

®T 

If  P>R  8t  =  Cos" HR/P) 

If  PiR  and  p  (of  the  navsat  closest  to  zenith) >1172; 
By  =  109.5°-  sin" HP  sin  19.5°/R) 

If  P*R  and p  (of  the navsat closest  tozenith)iir/2; 
0T  =  70.5° -sin" HP  sin  19.5°/R) 

DELONE 

If  Sin2  0T  >  C.  8j  ^Sin'Hc/SinB,) 

If  Sin2  0T  i  C,  8|  =0T 

DELTWO 

S2 

If  Sin2  8j  >  C,  82  =  8j 

If  Sin2  0T  i  C,  82=  Cos"Hl-2C)-8T 

NJL 

Number  of  navsats 

ISCMP 

Navsat  selection  technique  parameter  ( ISCMP =0 
for  zenith;  ISCMP=1  for  all  satellites  taken  four 
at  a  time) 

P(N,K) 

Array  containing  the  orbital  elements  for  the  satellites 
(where  N  is  the  identification  number  of  the  navsat  and 
k= 1  to  5  is  the  index  on  the  first  five  orbital  elements) 

PER(N) 

Orbital  period  of  navsat  or  user  satellite 

TABLE  XVI  (Continued) 


FORTRAN  EQUATIONS  EXPLANATION 


CK  C  Numerical  value  representing  the  fraction  of  the  area 

of  a  sphere  with  radius  equal  to  the  orbital  radius  of 
navsats;  used  in  calculating  width  of  band  where  best 
satellites  will  be  sought  ( must  be  set  equal  to  one  when 
navsats  are  in  elliptical  orbits) 


LATDEG 


ELEVAT 

LATIC.LONIC 

INC.ITF 

R,RMX 

UPV 

STN 

SE 

UTS  • 

USUV 

NSTO 

NSPL(L) 

CLU) 

Z 

YE 

XN 


Latitude  step  size  ( 5  or  10  degrees) 

Masking  angle  of  users 

Latitude  and  longitude  increments  to  be  used 

Time  Increment,  in  minutes,  at  which  calculations  8re  desired 
and  the  total  number  of  increments,  plus  one,  desired 

Vector  from  the  center  of  the  earth  to  a  satellite 

Vector  from  the  center  of  the  earth  to  a  user 

Vector  from  a  navsat  to  a  user 

Elevation  ( masking)  angle  of  a  satellite 

Vector  from  a  user  to  a  navsat  ( -STN) 

Unit  vector  (of  UTS) 

Total  number  of  navsats  in  view 

Total  number  of  navsats  in  view  at  each  latitude 

Total  number  of  latitudes  when  there  are  four  or  more  navsats  in 
view 

Vector  in  the  polar  direction  (origin  at  earth  center) 

Vector  in  the  eastward  direction  (origin  at  the  user) 

Vector  in  the  northward  direction  (origin  at  the  user) 


TABLE  XVII  (Continued) 


FORTRAN  EXPLANATION 


Q(  NSGD(  N),l )  Direction  Cosines  (where  NS0D(  N)  contains  the  “identification 

numbers"  of  the  four  selected  navsats  and  I  =  1  to  4) 

SIGT(N)  An  array  containing  the  DOP  parameters 

CDOP(L,K,IDOP)  Storage  for  DOPs  at  each  time  step;  L=  latitude  index, 

K=  longitude  index,  I  DOP  =  DOPs  index 

PIB(  IX)  Elevation  distribution:  the  probability  that  the  satellites  in  view 

will  have  specified  elevation  angles;  IX=  elevation  angle  index 

CAGX(  LA.LC)  Latitude  elevation  distribution:  the  probability  that  any  navsat  in 

view  will  have  an  elevation  angle  greater  than  or  equal  to  those 
specified;  LA=lat1tude  index,  LC=  elevation  range  index 

0LEB(  1C)  Accumulative  elevation  distribution:  the  probability  that  the 

elevation  angle  to  a  navsat  is  greater  than  or  equal  to  those  listed; 
IC=  elevation  range  index 

QSR(  IQSR)  Range  into  which  the  printed  variable  falls;  IQSR  *  1 ,36  in  steps 

of  .2 

SKEGX(  LK  ,IQSR  ,JDOP )  DOP  parameters  for  overall  global  performance;  LK=  latitude 

index,  IQSR  (see  above),  JD0P=  DOPs  index 

MAX(  IL.IK)  Maximum  number  of  navsats  seen  at  the  intersections  of  latitudes 

and  longitudes;  IL=latitude  index,  IK=  longitude  index 

MIN(  IL.IK)  Minimum  number  of  navsats  seen  at  the  intersections  of  latitudes 

and  longitudes;  IL=latitude  index,  IK=  longitude  index 


FORTRAN 


EXPLANATION 


OBLAT(IL,N) 

Probability  (in  percent)  of  seeing  exactly  N  navsats;  IL  = 
latitude  index,  N=number  of  navsats 

ACLAT(IL,N) 

Probability  (in  percent)  of  seeing  N  or  more  navsats;  IL  = 
latitude  index,  N=number  of  navsats 

OBDIS(N) 

On  a  global  basis,  the  probability  ( in  percent)  that  exactly 

N  navsats  will  be  seen 

ACTOT(N) 

On  a  global  basis,  the  probability  ( in  percent)  that  N  or  more 
navsats  will  be  seen 

NDEL 

Number  of  satellites  to  be  deleted  from  the  constellation 

JDEL(I) 

An  array  containing  the  "identification  numbers"  of  deleted 
satellites 

Explanation  of  Subroutines 

Table  XVIII  lists  the  individual  subroutines  utilized  in  the  program  and  provides  a  brief 
explanation  of  the  purpose  of  each.  A  more  detailed  explanation  of  each  subroutine  can  be  found  in 


Reference  1 7. 


TABLE  XVIII 


Computer  Program  Subroutines  ( 1 7:23-2?) 


SUBROUTINE  NAME 

DESCRIPTION 

ORBINI(N) 

Initializes  the  orbital  elements  for  the  navsats  and  user  satellite 
from  the  input  data 

ORB  1  T(  1  ,T  ,PER  ,R  ,VEL  ,AC) 

Iterates  for  the  eccentric  anomaly  and  computes  the  true  anomaly 
for  each  time  step 

TRMATX(TR.I) 

Calculates  the  3  x  3  coordinate  transformation  matrix,  TR 

MATMUL(T,Y,0) 

Performs  coordinate  transformation  of  a  vector  by  matrix 
multiplication 

POI  NT  ( ALO  ,ALA,T  1 M  ,VEC) 

Calculates  the  vector  from  the  center  of  the  earth  to  the  user  at 
a  specific  latitude  and  longitude 

TAAT(  MAX  ,MXX  .MATRIX) 

Sets  up  the  sequence  of  navsats  which  are  to  be  examined,  using 
either  the  one  above  or  below  the  user  as  one  of  the  four  in  each 
calculation  of  the  tetrahedron  volume 

VOLUME(  UVEC.IDSAT  ,VOL ) 

Calculates  the  volume  of  the  tetrahedron  formed  by  the  set  of 
four  satellites  used 

COVNAVf  G  ,1 D  ,NAT  ,SIG) 

Computes  the  six  dilution  of  precision  values  ( DOPs) 

BLOCK  DATA 

Contains  various  parameter  values  used  in  the  program 

VECT0R(Y1 ,1,V2,Y3) 

Performs  vector  additions,  subtractions,  and  cross  products 

DOT(VI  ,V2) 

Function  which  calculates  dot  product  of  two  vectors 

UNIYEC(V.UY) 

Calculates  unit  vectors 

TALL(  MAX  ,MXX  .MATRIX) 

Sets  up  the  sequence  of  navsats  which  are  to  be  examined,  using 
all  satellites  taken  four  at  a  time 

SATDEL(P,PER,NJ,ND,JDEL)  Deletes  specified  satellites  from  thenominal  constellation 


Appendix  C 

Computer  Program  Listing 


1. 

MAIN  PROGRAM 

[77-90] 

2. 

ORBINI(I) 

[91] 

3. 

ORBIT  (l,T, PER, R,VEL,AC) 

[91-92] 

4. 

TRMATX  (TR,I) 

[92-93] 

5. 

MATMUL  (T.V.O) 

[93] 

6. 

POINT  (ALOALA.TIM.VEC) 

[93] 

7. 

TAAT  ( MAX ,MXX  .MATRIX) 

[93] 

8. 

VOLUME  (UVEC.IOSAT.VOL) 

[93-94] 

9. 

COVNAV  (0.ID.NAT.SIG) 

[94] 

10. 

BLOCK  DATA 

[95] 

11. 

VECTOR  (V!,I,V2,V3) 

[95] 

12. 

DOT  (VI  ,V2) 

[95 

13. 

UNIVEC(V.UV) 

[95] 

14. 

TALL  (MAX, MXX.MATRIX) 

[95-96] 

15. 

SATDEL  ( P  ,PER  ,NJL  ,NDEL  ,JDE  L ) 

[96] 

’6 


THIS  PflOCBflH  IS  fl  MODIFICATION  OF  A 
THE  RAND  CORPORATION  ON  THE  GEOMETRIC 
RANGING  NAUIGRTION  SATELLITE  SVSTENS, 


SSES3E 


TIE  U8AF.  ** 


COHIOH/ORB I  S/P(37, 25 )/CON/C<  18  > 

DIMENSION  AA(3>,  8(37,4),  I  OS  (37),  ISIC(4),KXX<  10626,4),NS0D(37),R(3) 
1,RF(3),AnX(37,3>,SIGT(0>,U(3>,UPV(3>,USUU(37,3>,UTS(3>,UE(3),Z(3>, 
2  VE(3),XN(3>,  I  SAME  (37),  I  INT(37),UEL<37),  I0UT(3?),STN(3> 

DINENSION  QSR(36),CAG(  19, 18>,GLEB<  18),PIB<  18),ACLAT(  19,36), 

1  ACT0T(3C),0BLAT(  19,36>,060IS(36),CL(  19>,HIN(  19,36  ),NAX<  19,36), 

2  CAGX<  19,  18),NSPL(  19>,SKEG<  19,36,6),SKEGX<  19,36,6),GL0S<6,36), 

3  COGP<  19,36,6),  ISRPP(32),PER<36) 

REAL  LAT,  LONG,  LATDEG,LADG,  LAD 
DINENSION  JOEL (10) 

KTB-48 
I PRINT-0 
Z<1  >0. 

2(2  M). 

Z(3)-1.E+10 

rOEL-O 


LOOt;  USER  ON  SATELLITE 

L0O2;  USER  ON  GROUND  AT  SPECIFIED  LATITUDE  AND  LONGITUDE 
L0O3;  GLOBAL  CALCULATIONS 

READ  (5,*)  LOC,NJL,ISCMP,NOEL 
IF(NDEL.EQ.O)  GO  TO  50 
READ  (5,*)  (JDEL(N),N-1,NDEL) 

SO  GO  TO  (100,110, 120), LOC 

100  00  105  K-1,5 

READ  (5,*)  (P(N,K),OI,NJL) 

105  CONTINUE 


REflO  <3,*)  <PER<N),N-I,NJL) 

IF<M3EL.HE.0>  CflLL  SflTDEL  <P, PER, NJL.NDEL, JOEL) 
NJN-HJL+1 

READ  <5,*>  <P<NJN,K),K-1,5) 

REflO  <3,*>  PER<NJI1>,AIN,CK 
REflD  <5,*>  INC, ITF 
GO  TO  155 
C 

110  DO  115  K-1,5 

REflO  «,♦>  <P<N,K),N-1,NJL) 

115  CONTINUE 

REflD  <5,*>  <PER<N),N-1,NJL> 

IF  (WEL.NE.O)  CALL  SflTDEL  <P, PER, NJL, NOEL, JOEL! 
REflD  (5,*)  flTL,ONGL,ELEURT 
REflD  <5,*>  INC, ITF 
00  TO  153 
C 

120  DO  125  K-1,5 

REflD  «,♦>  <P<H,K),H-1,NJL> 

125  CONTINUE 

REflD  <5,*>  <PER<H>,N-1,HJL> 

IF(NDEL.NE.O)  CflLL  SflTDEL  <P, PER, NJL, NOEL, JOEL) 

REflD  <5,*>  LflTDEG,ELEURT 

REflD  (3,*)  LflTIC,LONIC, INC, ITF, IPFflEQ, ITIHE 

00  135  Hfl-1,10 

CL<NRM). 

NSPL<HRM) 

DO  130  IS- 1,36 

HIN<HR,NB)-30 

NRX<NR,NBM) 

DO  130  HC-1, 6 
SKEGX(NR,NB,NCM). 

coop<HR,re,nc)-o. 

130  SKEO<NR,ie,HC)H). 

DO  135  tf)-1, 18 
CflOX<Nft,NDW). 

135  Cfl0<HR,re>-0. 

DO  140  HOI, 36 

oo  140  re-i,6 
140  GL0S<re,NR)-0. 

DO  130  Hfl-1,10 
RCT0T<HR)-0. 

oboisovd-o. 

00  145  ND-1, 36 
RCLflT<NR,re)-0. 

143  O0LflT<HA,HDM). 

150  0L£B(Hfl)-O. 

C 

155  CONTINUE 
PRINT  105 
PRINT  100 
00  160  IP- 1, NJL 

PRINT  213,  IP,<P<IP,IU),IU-1,5),PER<IP) 

160  CONTINUE 

IF<LOC.EQ. 1)  PRINT  200 
IP-HJL+1 


IF<LX.EQ.  1 )  PRINT  215,  IP,<P<IP,  IU>,  IU-1,5),PER(IP> 
IFCL0C.EQ.2)  PRINT  210,  ATL,OHGL,ELEURT 


IFCL0C.EQ.3)  PRINT  205,  ELEUAT,LATOEG,LATIC,LONIC, I T I ME 
ITT»<ITF-1)*INC 
PRINT  220,  ITT, INC 
IF<L0C.EQ. 1)  PRINT  225,  AIN,CK 
IF<L0C.EQ.  1.RND.  ISCNP.EQ.O)  PRINT  100 
IFGLOC.NE. 1.AND. ISCHP.EQ.O)  PRINT  105 
IF<I8CNP.EQ. 1>  PRINT  170 
C 

C  TIE  FOLLOUING  FORHRTS  HflUE  TO  DO  WITH  INPUT 
C 

105  FGRHATdHO,  10X,  ’THE  SRTELLITE  HOST  NERRLV  OUERHERO  IS  USED  RS  ONE 
10F  THE  FOUR'/,  1 1X,  *  IN  ALL  CALCULATIONS  OF  THE  UOLUTC  OF  THE  TETRflH 
2EDR0N*  > 

100  FORNRTC 1H0, 10X,  'THE  SRTELLITE  HOST  NERRLV  ABOUE  OR  BELOU  IS  USED  R 
IS  ONE  OF  THE  FOUR'/,  1 1X, '  IN  ALL  CALCULATIONS  OF  THE  UOLUHE  OF  TIE 
2TETAAHEDA0N*  > 

170  FORHRT<  1H0, 10X,  'ALL  SATELLITES,  TAKEN  FOUR  AT  A  TINE,  ARE  USED  IN' 
1,/,11X,*1HE  CALCULATIONS  OF  TIE  UOLUHE  OF  TIE  TETRAHEDRON') 

175  FORNRT<  1015) 

180  FORHAT(  12F0.0) 

105  FORHAT<7F10.0> 

100  FORHAT< 1H0,22X, '  ECC* ,5X, 'ARGP' ,4X, 'RASC* ,4X, ' INC' ,5X, 'ANON' ,4X, 
1‘PER',//) 

195  FORMAT < 1H1, //32X, '0R8ITRL  ELEMENTS') 

200  FOAMRTC  1H0,10(X,  'USER  SATELLITE  ORBITAL  ELEIENTS’/) 

205  FORHAT(  1H0, 10X,  ‘GLOBAL  DISTRIBUTION  CALCULATIONS' ,  //, 

1  1 1X, 'HRSKING  ANGLE  «  ',F5.2,'  DEGREES',/, 

2  1 1X,  'LATITUDE  STEP  *  ',F6.2,'  DEGREES',/, 

3  11X, ‘LATITUDE  INCREMENT  ■  ',13,/, 

4  1 1X, 'LONGITUDE  INCREMENT  -  ',13,/, 

5  1 1X, ‘DILUTION  OF  PRECISION  PARAMETERS  PRINTED  AT  TIME  INCREMENT  0 
OF  ',15) 

210  FORMAT <  1H0, 10X,  'USER  LOCATION  ON  EARTH' , /I IX, 

1  ‘LATITUDE  -  *,F5.2, '  DEGREES' /,  1 1X,  'LONGITUDE  -  ',F5.2, 

2  '  DEGREES',/,  11X, 'HRSKING  ANGLE  «  ',F6.2,'  DEGREES') 

215  FORHAT< 1H  , 15X, 13, 1X,F0.3,5F0.2) 

220  FORHRT<  1H0, 10K,  'TOTAL  TINE<HIN)  ■  *,  15,  /,  1 1X, 

1  'TINE  INCREHENT(HIN)  -  ',13) 

225  FORHRT< 1H  , 10X, ‘BEAMUIDTH  ANGLE (DEG)  -  ' ,F6.2, /, 1 IX, 

1  'FRACTION  OF  NAUSRT  SPHERICAL  AREA  =  ',F5.3> 

C 

C  SET  UP  ORBITAL  ELEMENTS 
C 

00  230  N-1,NJL 

P<H,22 >-22000. *<PER<N>/24.  W2./3.  > 

CALL  ORBINI  (N) 

230  CONTINUE 

IF<L0C.EQ.2.0R.L0C.EQ.3>  GO  TO  235 
P<IP,22«2000.*<PER<IP>/24.  W2./3.  > 

CALL  ORBINI  <IP> 

235  CONTINUE 
INCAM) 

N8TOO 

HRXNSS-0 

DO  025  IT-1, ITF 

TIO^L0ATCINCA)/144O. 

ITOUT-< IT-1 )*INC 
INCA-INCAHNC 


IF(L0C.EQ.2.0R.L0C.EQ.3>  ftHUL 

oo  240  m\tm 

CALL  ORBIT  CN,TJD,PER,AF,UE,RA> 
00  240  10*1,3 
RNXCN,IU)-flFCIU> 

240  CONTIHUE 

IFCLQC.EQ.2)  GO  TO  335 
IFCL0C.EQ.3)  GO  TO  330 


DO  245  10*1,3 
UPVCIUMWXCIP,IU> 

245  CONTINUE 

IK-0  _ 

□cm  in  4.  oooooo 
00  205  N*1,NJL 
DO  250  10-1,3 
RCIU)-RHXCN,IU> 

250  CONTINUE 

CALL  OECTOR  <R,2,UP0,0T3) 

AR-SQRTCDOTCR,R» 

RP*SQRT<DOT<UPO,UPO)) 

B-RCOSC2.O920144E+O7/AR) 

IFCRP.LE.2.0925144E+0?)  PRINT  251 

251  F0RHRTC1H  ,/,3X,  TERMINATION  OF  RUN,  ALTITUDE  APPROACHING 
IFCAP.LE.2.0920144E+07)  STOP 

B-RC06C2 . 0020 144E+07/AP  ) 

PHI*R*B 

OOTRP-OOTCR.UPU) 

TVETN-RCOS<DOTRP/CRR*flP>) 

IFCTHETN.GE.PHI)  GO  TO  209 
RLHSGRT  (DOT  CUTS,  UTS  )  ) 

DOmjP-OOT<UTS,UPO> 

BETArMCOSC-OOTUP/CRU*AP )  ) 

AIUID*AIN/C<2> 

BHIOTINCOMTHETmflETAN) 

IFCRIUlD.LT. BUIDTH)  GO  TO  205 

IK-IK+1 

ISRUECIKMf 

BETRSO-GETRN 

IFCBETAN.GT.C<5>>  BETATKO >-BETAN 
255  IFCBETflN.LT. BETATtO  GO  TO  200 
GO  TO  209 
200  I8ATN0-N 
BCmtNBETAN 
BETRSHBETR80 
209  CONTINUE 

IFCIK.LT. 4>  GO  TO  022 
00  285  1-1, IK 

IFCISRUECD.EQ.ISRTNO)  GO  TO  270 
GO  TO  285 
270  INOX- I 

00  275  KI-1, IHDX 
KII-KI+1 

I INTCKI I  XSRUECKI  > 

275  CONTINUE 
INBX-INDX+1 


oooo 


I  INT<KI  >ISRUE<KI  > 

280  CONTINUE 
00  TO  280 

289  CONTINUE 

290  IINT<1)-ISATN0 
IFdSCNP.EQ.  1)  00  TO  339 

TIC  NUTfiERS  OF  TIC  SATELLITES  WHICH  FIT  CRITERIA  FOR  USE,  PLUS  THE 
OTC  NEAREST  TO  OUERCflD  HRUE  BEEN  CALCULATED 

IOS< 1  >■! INT<1> 

NSS>1 

DO  295  10-1,3 
R<iu)=flnxcMmxi>,  iu> 

299  CONTINUE 

CALL  UECTOR  <R,2,UPU,UTS> 

CALL  UNIUEC  <UTS,U> 

DO  300  10*1,3 
USUU<I  INTO),  IVMKIU) 

300  CONTINUE 

DO  329  N-2,  IK 
DO  309  10*1.3 
ROUVANXCI  INT<N),  IU) 

305  CONTINUE 

CALL  UECTOR  <R,2,UPV,UTS> 

CALL  UNIUEC  <UTS,U> 

DO  310  10-1,3 
USUKI  INT<N),  IUMKIU) 

310  CONTINUE 

AR-SqRT(DOT(R,R» 

AP*SQRT(DOT(UPU,UPU)> 

IF<AP.OT.AR>  THETT-RCOS<AR/AP> 

IF(AP.LE.AR. AND. BETAS. GT. C(9)>  THETT- 109. 9/C<2  MIS  I  N<AP*.  33380660' 
1AR) 

IF<AP.LE.AR. AND. BETAS. LE.C<9»  THETT»70.9/C<2MISIH<AP*.33380060/B 
1R> 

DOTRP-OOT<R,UPO> 

THETTMC08<D0TRP/<AR*RP  >  > 

SINTT-SIN(THETT) 

SINTSQ-SIHTT**2 
IF<SINTSQ.GT.CK>  GO  TO  319 
OELONE=THETT 
ANt^l  .-<CK/.3) 

DELTUO-ACOS(ANG J-THETT 
00  TO  320 

319  CONTINUE 
DELONE-RSINCCK/SINTT) 

DELTMMELONE 

320  IF<THETN.LT.(THETT-OELONE>)  00  TO  329 
IFCTHETN.GT. <THETT+OELTUO)>  GO  TO  329 
H8S-N88+1 

I08<N88>I  INT<N) 

329  CONTINUE 

IF<NSS.LE.4)  GO  TO  629 


':*■  «v'-,  »  t  -  »■. 


L  w  n  *  V  ■  -  •  w  -  •  -  V 


C  END  OF  USER  ON  SATELLITE 
C 

C  CALCULATIONS  FOR  USER  ON  EARTH  AT  SPECIFIED  LflT  AND  LONG,  AND 
C  GLOBAL  DISTRI0UTIUON8  FOLLOW 
330  CONTINUE 
LKL-36 
LLL-10 

339  CONTINUE 

IFCL0C.EQ.2.0R.LQC.EQ.1)  LKL-1 

IF<LOC.EQ.2,OR.LOC.EQ. 1 >  LLL-1 

DO  900  K«1,LKL,L0NIC 

DO  900  L»1,LLL,LATIC 

IF<LOC.EQ. 1.AND. ISCNP.EQ.O)  00  TO  401 

IFCLOC.EQ. 1.AND. ISCNP.EQ. 1)  GO  TO  379 

LQNG-FLDRTCK-1>*10. 

LAT-90. -FLOAT CL- 1  )*LATI3EG 
IFCLOC.EQ. 2)  LONG-ONGL 
IFCLOC.EQ. 2)  LAT-ATL 
CALL  POINT  CLONG,LAT,TID,UPU> 

NSS-0 

CLCLXLCLH1. 

00  390  N"1,NJL 
00  340  IU-1,3 
RCIUWWXCNJU) 

340  CONTINUE 

CALL  HECTOR  CUPV,2,A,STH> 

SE— DOTCSTN,  UPU  )/SQATCDOTCSTN,  STH  >*OOTCUPV,  UPV  » 
IFCAB8CSE).GE. .0009000)  SE«8IGNC1.,SE> 

ELaA8INCSE)*CC2> 

IFCEL.LT. ELEURT)  GO  TO  390 
NSS-NSS+1 
I INTCNSSM1 
UELCNJ-EL 

CALL  UECTOR  CR,2,UPV,UTS> 

CALL  UNIUEC  CUTS.U) 

DO  349  10*1,3 
U8UUCN,  IVMJCIV) 

349  CONTINUE 
390  CONTINUE 
399  CONTINUE 

NSTO-NSTOtNSS 
NSPLCL  MISPLCL  HNSS 
NSP-NSS+1 

HRXN8S-HRX0CHRXNSS,NSS ) 

HIGH-O. 

DO  369  NUN-1,NSS 
IX-I INTCNUH) 

AEfMJELCIX) 

IFCAEA.GT.HIGH)  GO  TO  360 
IFCAEA.LE.HIOH)  GO  TO  369 
360  HIGH-AEA 
NN-NUN 

309  CONTINUE 

00  370  NU-1,NSS 
IFCNU.EQ.NN)  NX-1 
IFCNU.LT. NN)  NX-NU+1 
IFCNU.GT.NN)  NX-NU 
I  OS  <  NX)- 1 INTCNU) 
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oooouooo 


370  CONTINUE 
GO  TO  400 

lOS(NSS)  HRS  THE  SATELLITE  NEWEST  TO  OUERHEAD,  THEN  ALL  OTH 
WICH  FIT  CRITERIA 

NEXT  CALCULATIONS  ARE  CONCERTED  WITH  FINDING  THE  COMB I NATION 
FOUR  SATELLITES  UHICH  HAUE  THE  GREATEST  UALUE  OF  THE  UOLUHE 
TETRAHEDRON  FORMED  BV  THEN 

379  CONTINUE 

DO  300  IKN»1, IK 


R<IU>ANX<I  INTtSKM),  IU> 
300  CONTINUE 

CALL  UECTOR  <R,2,UPU,UTS> 
CALL  UN I  ICC  <UTS,U> 

DO  385  IU»1,3 
USUUd  INTdKH),  IUMKIU) 
385  CONTINUE 
300  CONTINUE 

DO  305  IJK-1JK 
305  IOS(IJK>Bl INTdJO 
00  TO  401 

400  IFCNSS.LT.4)  GO  TO  470 


IF(ISCMP.EQ.I)  GO  TO  405 

K0T-N8S 

NSS-HSS-1 

CALL  TART  <NSS,KCON,KXX> 
ISIC(1>I0S<I> 

LPN-2 
GO  TO  410 
KOT-IK 

IF<LOC.£Q.2.0R.LOC.EQ.3> 
CALL  TALL  <KOT.KGCN.KXX> 


00 

DO  415  LPQ-LPN,4 
NUX*400C<H,LPQ> 
ISIC(LPQVIOS<NUX) 

415  CONTINUE 

CALL  UOLUNE  <USUU, ISIC,UOLUM> 
IF<UOLUN.GT.BOX)  GO  TO  420 
IF(UOLUN.LE.BOX)  GO  TO  430 


DO  425  Ht>1, 4 
NSGO(HC)-ISICCHC) 
425  CONTINUE 
430  CONTINUE 


DO  435  101,3 
R(IU>AMX(NSGO<N),  IU> 

435  CONTINUE 

CALL  UECTOR  (R,2,UPV,UTS> 
CALL  UECTOR  <Z,3,UPU,VE> 
CALL  UECTOR  (UPU,3,VE,XN> 
AX-8QRT<00T<XN,  70) 


AV-SQRT<DOT<VE,VE» 

AP-SQRT<XT<UPV, UPU  >  > 

RLT-RP/6076. 116-3444. 

AS-SQRT<XT<UTS,UTS>> 

Q<NSX<N  >,  1  >»OOT<XN,  UTS  >/<AX*RS  > 

G<NSX<N  ),  2  >-XT<  VE,  UTS  >/<AV*RS  ) 
G<NSX<N>,3«0T<UPU,UTSV<AP*flS> 

G<NSX<H>,4>1. 

440  CONTINUE 

CALL  COUNflW  <G,NSX,4,SIGT> 

IF(LX.EQ.3>  X  TO  480 
IF(L0C.EQ.2>  ALT-O. 

X  443  1 1-1,4 
IOUT<IO<NSGD<I  I  > 

445  CONTINUE 
IL-1 
LK-4 

DO  455  NN-ljKOT 
ITST=NSGO< IL) 

IFOTST.EQ.  IOS<NN>>  X  TO  450 
LK-LK+1 

IOUTOJO=IOS<NN> 

X  TO  455 
450  IL*IL+1 

IF< IL.EQ-5)  X  TO  460 
455  CONTINUE 
460  KOS-HN+1 

X  405  LK»KOS,KOT 
IOUT<UO-IOS<UO 
465  CONTINUE 

ITOUT-< IT-1 >*INC 
KTR=KTTH2 

IF<K0T.GT. 18)  KTB-KTR+1 
IF<KTR.GE.40)  X  TO  470 
X  TO  475 
470  PRINT  525 
KTB-0 

475  PRINT  530,  ITOUT,flLT,<SIGT<KP),KP»1,6),<IOUT<KP),KP»1,KOT> 
C 

C  END  OF  CALCULATION  FOR  SINGLE  USER  ON  ERRTH 
C 

X  TO  625 
C 

C  FOLLOWING  CUMULATIONS  FOR  GLOBAL  DISTRIBUTION 
C 

479  KOT-NSS 

SIGT<1)-10X. 

SIGT<2)-10X. 

SIGT<3>10X. 

SIGT<4>-10X. 

SIGT<3)-10X. 

SIGT<6)»10X. 

4X  X  485  I00P«1,0 

485  COOP<L,K,IOOP)-SIGT<IOOP> 

X  490  NS- 1,6 

KR-fVlXD<1,niN0<INT<SIGT<NS>*5.+1.>,36» 

4X  SKE0<L,KR,NS)-SKEG<L,KA,NS>t1. 

X  495  NR-1,K0T 


84 


I-IOS(NA) 

ELT*UEL< I ) 

1(0=41 1  N0(INT(ELT/5.  HI,  18) 

485  CflOCL,  KO  HCttKL,  KO  )+ 1 . 

NRX(L,  K)=NAXO(NAX(L,  K  ),  KOT  ) 
t1IHa,K)-t1INO((1IH(L,K),KOT) 

0BLRT(L,  MSP  >0BLfiT<L,  NSP  H 1 . 

500  CONTINUE 

IF(L0C.EQ.3.RND. IPFREQ.EQ.O)  00  TO  625 
I  NO  INTERNED I RTE  PRINT 

IFdPRINT.EQ.O.RND.  IT.EQ.  1)  GO  TO  505 

PRINT  FIRST  TINE  STEP 

IPRINT=IPRINT+1 

IF< IPflINT.EQ. ITINE)  GO  TO  505 

PRINT  ERCH  TINE  STEP  REQUESTED 

GO  TO  625 
505  IPRINT-0 

DO  515  ID0P*1,6 
IF(IDOP.EQ.I)  PRINT  535,  (TOUT 
IF( ID0P.EQ.2)  PRINT  540,  I TOUT 
IF<ID0P.EQ.3)  PRINT  545,  ITOUT 
IF< ID0P.EQ.4)  PRINT  550,  ITOUT 
IF(IOOP.EQ.S)  PRINT  555,  ITOUT 
IF( ID0P.EQ.6)  PRINT  560,  ITOUT 
IF<LRTDEG.EQ. 10. >  PRINT  565 
IF(LRTDEG.EQ.5. )  PRINT  570 
ICT=-1G 

DO  510  101,36 
ICT-ICT+10 

IF< ICT.LE.90.0R. ICT.GE. 190)  PRINT  575,  ICT, (CDOP( IK, 1C, IDOP>, IK 
1*1,19) 

IFdCT.EQ.  100)  PRINT  580,  (COOPdK,  1C,  IDOP),  IK*  1, 19) 

IFdCT.EQ.  110)  PRINT  585,  <COOP< IK,  1C,  IDOP),  IK-1, 19) 

IFdCT.EQ.  120)  PRINT  590,  <CDOP<IK,  1C,  IDOP),  IK*1, 19) 

IFdCT.EQ.  130)  PRINT  595,  <COOPdK,  1C,  IDOP),  IK=1, 19) 

IFdCT.EQ.  140)  PRINT  600,  (COOPdK,  1C,  IDOP),  IK=1, 19) 

IFdCT.EQ.  ISO)  PRINT  605,  (COOPdK,  1C,  IDOP),  IK*1, 19) 

IFdCT.EQ.  160)  PRINT  610,  <COOP<  IK,  1C,  IDOP),  IK-1, 19) 

IFdCT.EQ.  170)  PRINT  615,  (COOPdK,  1C,  IDOP),  IK*1, 19) 

IFdCT.EQ.  180)  PRINT  620,  (COOPdK,  1C,  IDOP),  IK-1, 19) 

510  CONTINUE 
515  CONTINUE 

DO  520  ia«1,19 
DO  520  ICK-1,36 
DO  520  I CD-I, 6 
520  COOPdCL,  ICX,  ICD)-0. 

THE  FOLLOWING  FORMATS  HRUE  TO  DO  UITH  R  SPECIFIED  TINE  STEP 
REQUEST  FOR  PRINTING 


525  FORMAT ( 1H1,/,  IX,  'TIIC(ffl)' ,  IX,  'flLT(ltl)'  ,4X,  UOOP' , 

1  5X,  ’HD0P‘,5X,  *H)0P',5X,  "TDOP*  ,5X,  'POOP* ,5X, 'GOOP' ,4X, 

2  'SATELLITES  CHOSEN',//) 

530  FORHRT(  1H0, 17,F9.0,6(F9.3>,  IX,  19(  IX,  12),  /,72X,  17(  IX,  I2» 

531  FORHRTdHO,  I7,F9.0, '  ONLV  ',14,'  SATELLITES  RRE  UISIBUE.') 

535  FOFMRTdHI,//,  10X,  'TINE  »  ',  I6,20X, 'UOOP  -  ALTITUDE'  ) 

540  FORNRTdHI,//,  10X,  'TINE  -  *,  I6,20X,  HDOP  -  POSITION  ERAOR  IN  HORIZ 


2/3 


RD-A151  692 
UNCLASSIFIED 


GLOBAL  POSITIONING  SVSTEN  -  A  MODIFICATION  TO  THE 
BASELINE  SATELLITE  CONS..  <U)  AIR  FORCE  INST  OF  TECH 
URIGHT-PATTERSON  RFB  OH  SCHOOL  OF  ENG I.  D  U  THOMIN 
DEC  84  AFIT/GA/ENG/84D-4  F/G  17/7  NL 


Hi  \  iij 

pH 

g^iiYJ 

pH 

IfriYj 

10HTRL  PLRME* ) 

545  F0lt«T<1H1,//,1OX,*TI1E  -  *, 

550  F0FmT<1H1,//,1OX,*TIHE  ■  *, 

555  FORrfrraHi,//,  iox,*time  «  \ 
non  ERROR*  > 

560  F0flmT<1H1,//,1OX,TirE  -  ', 
ion  ERROR*  > 

565  FOfinflrmHO,52X/LflTITlJOE,,//)«<  1  90*, 3X 
1  3X, *  50*, 3X, *  40’, 3X, *  30*,3X,*  20*, 3X 
3X,  *-30*, 3X,  *-40*,3X,  *-50*, 3X 
/> 

52X, ‘LflTmJ0E*,//,9X, *  90*, 3X 
3X, *  65*,3X, *  60* ,3X, *  55*, 3X 
3X. *  30‘ ,3X, *  25* ,3X, *  20*,3X 


575  FonmT<tH 
560  FOnmidH 
565  FOmmiH 
590  FORHRT< 1H 
595  FanmT<iH 
600  FORHRK 1H 
605  FOnmTdH 
610  FommaH 
615  FamraH 
620  FORNRfK  1H 


,2X, 13, IX, 19F6.2) 

, *L*, IX, *  100*, IX, 19F6.2) 
,*0*,1X,*110*,1X, 19F6.2) 
, *M*, IX, *  120*, IX, 19F6.2) 
,‘G*,1X,*1W,1X,  19F6.2) 
, *  I *, IX, *  140*, IX, 19F6.2) 
, *T‘, IX, *  150*, IX, 19F6.2) 
,  'll*,  IX,  *  160*,  IX,  19F5.2) 
,  'O',  IX,  *  170*,  IX,  19F6.2) 
, *E*, IX, *  180', IX, 19F6.2) 


622  IT0UT«<IT-1)*irC 

RLT-flP/0070. 1 16-3444. 
KTTHCTFH2 

IF  <KTR.GE.40)  00  TO  623 


623  PRINT  525 
KTR-O 

624  PRINT  531,  ITOUT,RLT,  IK 

625  CONTINUE 


DIHENSIQNRL  POSIT 
DINENSIONRL  P08ITI 


3X, ' 

70* 

,3X 

3X,  * 

0* 

,3X 

3X,* 

-70* 

,3X 

3X, ' 

80* 

,3X 

3X, ' 

45* 

,3X 

3X,* 

10* 

,3X 

630  QSR<LI)-FL0RT<LI-1>*.2 
DO  640  IX-1, 18 
PIB<IXM >. 

DO  635  IY-1,19 

635  PIB<IXM>IB<IX>*CRQ<IV,IX> 

640  PIB<IXXPIB<IX)/FL0RT<N8T0))*100. 

DO  661  Lf*1,19 

IF<CL<LR>.EQ.O..OR.NSPL<LA>.EQ.O>  00  TO  I 
DO  690  ft-1,36 
DO  650  1-1,6 
9KE0X<Lfl,N,  I  M). 

DO  645  J-N,36 

645  8KE0X(LR,N,  I  >-SKEOX(LR,N,  I  >+SKEG<LA,J,  I  > 
650  SKEGX<LR,N,  I  >*SKEGX<LR,N,  I  >/CL<LR> 

DO  660  IX-1, 18 
CR8X<LA,LCM). 

00  695  LF<C,  18 

695  CRQX(LR,LCXfleX<LR,LC>4CR0<LR,LF) 

CR0X<LR,LXXR0X<LR,IX>/FL0RT<N8PL<LR» 


660  CONTINUE 

661  CONTINUE 


PRINT  780 

PRINT  785,  <PIB<  IX),  IX-1, 18) 

IFCLRT0EG.EQ.5. )  PRINT  795 
IFOATDEO.EQ.  10.  >  PRINT  800 
ICT— 5 

00  665  LC-1,18 
ICT-ICT+5 

PRINT  815,  ICT, <CAGXOfl,LC),LR-1, 19) 

665  CONTINUE 
CUS-O. 

00  680  l.n-1, 19 

IF<a<LH>.EQ.O..OR.NSPL(U1).EQ.O>  00  TO  680 
CIM»S«9O.-FL0ATttJ1-1  PlflTDEG  >*C<  1 » 

CUS-CUS+CU 
00  670  IN-1,6 
DO  670  11-1,36 

670  GLOB<  IN,  1 1 MLOSC  IN,  1 1  >+SKEGX<LH,  1 1 ,  IN)*CW 
00  675  LN-1,18 

675  GLEB<LN>0LEB<U1>*CflGX<Li1,LN)*CU 
680  CONTINUE 

00  685  LN-1, 18 
685  OLEB(LN  )-GLEB<LN  )  /CUS 
PRINT  805 

PRINT  810,  <GLEB<  1C),  101,18) 

00  690  IN-1,6 
00  690  11-1,36 

690  GLOSCIN,  1 1 MXJDSOH,  1 1  )/CUS 
DO  695  JDOP-1,6 
IF(JDOP.EQ.I)  PRINT  830 
IF<JD0P.EQ.2)  PRINT  835 
IF<JD0P.EQ.3)  PRINT  <M0 
IF<JD0P.EQ.4)  PRINT  845 
IFCJD0P.EQ.5)  PRINT  850 
IF<JD0P.EQ.6)  PRINT  855 
IFOflTDEG.EQ.  10. )  PRINT  820 
IFOflTDEG.EQ.  5. )  PRINT  825 
00  695  IQSR-1,36 

PRINT  775,  QSROQSR),  (SMEBX<LK,  IQSR,JD0P),LK-1, 19) 
695  CONTINUE 
PRINT  860 
00  700  IQSR-1,36 

PRINT  865,  QSR( IQSR),<OLOS(  IN,  IQSR),  101,6) 

700  CONTINUE 
PRINT  870 
LA0090. 

IFOATDEG.EQ.tO.)  LAO- 10. 

IFOflTDEO.EQ. 5.)  LFO-5. 

PRINT  880 
00  705  IL-1, 19 

PRINT  885,  LflDO,(NRX(IL,  IK),  IK-1,36) 

LAOOLADOHflD 
705  CONTINUE 
PRINT  875 


i 


LRDO-QO. 

00  719  IL-1, 19 
00  710  IX-1,36 

IF<HIH<IL,  IX).EQ.30>  HIN<IL,IX)-0 
710  CONTINUE 

PRINT  889,  U»G,<HIN<IL,  IK>,  IK-1, 36> 
LflOG-LflOGHJC 

719  CONTINUE 

DO  720  N-1,I*SSP0 
060IS<N)*0. 

DO  720  L-1,19 

IF<CL(L>.EQ.O.  .OR.NSPL(L).EQ.O)  00  TO  720 
06LRT(L,NXOBLRT<L,N>/CL(L>)*100. 

720  CONTINUE 

DO  730  L-1, 19 
00  730  N>1,NNSSF0 
RCLBT<L,N>0. 

DO  729  1H1,mSSP0 

729  RCU1T(L,N>4)CLRT<L,NH06LRT(L,H> 

730  CONTINUE 
CM). 

00  740  L-1, 19 

IF<CL(L).EQ.O. .OR.NSPL(L).EQ.O)  80  TO  740 
CfKOS(  <90 .  -FLOfiT  <L- 1  *LflTDEG  >*C<  1 » 
COCO+CR 

oo  739  n*i,msspo 

739  080IS<N>-aeDIS<N>*aeLfiT<L,N>N» 

740  CONTINUE 

DO  745  N-1,mSSP0 
749  080IS<NMM)IS<N)/C0 
DO  790  N-1.NNSSP0 
fiCTOT<HM>. 


rctot<n  wctot<m  >+oeo  i  s<n  > 

00  799  1-1,32 

ISflPP<l>«l-1 

CONTINUE 


NNMMSSPO 

IF  <mSSP0.6T.16>  HNS- 16 
PRINT  890,  <ISflPP<l  >,  1-1, 16) 

PRINT  900 
LRDG-90. 

00  769  IL-1,19 

PRINT  899,  LflD0,<08LflTCIL,H>,rH#>,mS) 

LflDG-LflDG+UC 

CONTINUE 

PRINT  909 

LRDO-90. 

DO  770  IL-1. 19 

PRINT  899,  LAOG,<RCLRT<IL,N>,N-NP,HNS> 

LflOMJBQHRD 

CONTINUE 

PRINT  910 

PRINT  790,  <060 1 S<N), IMF, HNS) 

PRINT  919 

PRINT  790,  <flCT0T<N>,fHf>,lt1S> 


IF(MNSSPO.LT. 16 >  GO  TO  920 
IF<MNS.GT. 16)  00  TO  920 

ms-mssro 

IP-17 

PRINT  890,  <I8RPP< I  ),  1-17,32) 

00  TO  760 
C 

C  TIC  FOLLOUINO  FORHRTS  HRUE  TO  00  UITH  fl  QLOBRL  SVSTEN 
C 

779  FOANATdH  ,F4. 1,3X,  19F6.3) 

780  FORNRTdHI, /////IX, 'ELEURTION  DISTRIBUTION  -  PROBRBILITV  THfiT  THE 
1SRTELLITE  IN  UIEU  HILL  HRUE  ELEURTION  ANGLES  R8  LISTED' ,//,48X, 

2  'ELEURTION  ANGLE*) 

789  F0RNRT(1H0,6X,  *  0-9  9-10  10-19  19-20  20-25  29-30  30-39  39-40  40 

1-49  49-90  90-99  99-60  60-69  69-70  70-79  79-80  80-89  85-90', //,2X, 

2  ‘PROS',  IX,  18F6. 1) 

790  FORHRT<  1H0,3X,  'PROS'  ,9X,  16F7.2) 

795  FORHAT<  1H0,  ////,  IX,  ‘LATI1U0E  ELEURTION  DISTRIBUTION' ,  /, 

1  IX, 'PROBABILITY  THAT  ANY  8RTELLITE  IN  UIEU  UILL  HRUE  AN  ELEUAT 
2IQN  ANGLE  GREATER  THAN  OR  EQUAL  TO  THOSE  LISTED', //,52X, 
3‘LRTITU0E',/,2X/ELEU‘,/,  IX, 'ANGLE',  IX, '  90  89  80  79 

47069809590494039302520 
5  15  10  5  0',/> 

800  FORHRTdHO,////,  IX, ‘LATITUDE  ELEURTION  DISTRIBUTION',/, 

1  IX,  ‘PROBABILITY  THAT  AW  SATELLITE  IN  UIEU  UILL  HRUE  AN  ELEURT 
2I0N  ANGLE  GREATER  THAN  OR  EQUAL  TO  THOSE  LISTED‘,//,92X, 
3'LRTITU0E',/,2X,'ELEU',/,  IX, 'ANGLE',  IX,'  90  80  70  60 

4  90  40  30  20  10  0  -10  -20  -30  -40  -50 


9  -60  -70  -80  -90',/) 

805  FOAHRT<  1H0, ////, IX, ‘ ACCUHULATHE  ELEURTION  DISTRIBUTION',/, 

1 1X, ‘PROBABILITY  THAT  THE  ELEURTION  ANGLE  IS  GREATER  THAN  OR  EQUAL 
2T0  THOSE  LISTED',  //,48X,  'ELEURTION  ANGLE') 

810  FORHRTdHO,  7X,'  0  5  10  15  202530354 

10  45  90  95  60  69  70  75  80  85',//,2X, 

2  'PROB' ,  IX,  18F6. 1 ) 

815  FORNRTdH  ,  13, 4X,  19F6.2) 

820  FORNAT<  1H0,52X,  'LATITUDE' ,  //,2X,  'NUN'  ,5X, 

1  '  90', 3X,’  80', 3X,'  70', 3X, '  60', 

2  3X, '  90', 3X,'  40’, 3X,'  30', 3X,'  20', 3X, '  10', 3X, '  0',3X,'-10', 

3  3X, '-20', 3X, '-30', 3X,  -40',3X, '-90', 3X,  -60', 3X, '-70', 3X, '-80', 

4  3X/-90',/) 

825  FORHRTdHO, 52X,  'LATITUDE', //,2X,  'NUN'  ,5X, 

W  ,3a,  ,3a,  W  ,3a,  73  , 

2  3X, '  70', 3X,'  69' ,3X, '  60', 3X,'  95', 3X, '  90’, 3X,'  45',3X, '  40', 

3  3X, '  39', 3X,'  30', 3X, '  29', 3X,'  20', 3X,'  15', 3X,'  10', 3X, '  9', 

4  3X, '  0',/) 

830  FORHATdHI,///,  IX, ‘DILUTION  OF  PRECISION  -  RCCUNULATIUE  LATITUDE  0 
1ISTRIBUTI0N',/,  IX, 'PROBABILITY  THAT  ALTITUDE  OOP  UILL  BE  GRERTER  T 
2HAH  HUNGER  LISTED',//) 

835  F0RHRT<1H1,///,  IX, ‘DILUTION  OF  PRECISION  -  RCCUNULATIUE  LATITUDE  0 
1I8TRIBUTI0N',/,  IX, ‘PROBABILITY  THAT  POSITION  OOP  IN  HORIZONTAL  PLA 
2HE  UILL  BE  GRERTER  THAN  NUMBER  LISTED',//) 

840  FORNRTdHI,///,  IX, ‘DILUTION  OF  PRECISION  -  RCCUNULATIUE  LATITUDE  0 
1ISTRIBUTI0N’,/,  IX, ‘PROBABILITY  THAT  LARGER  COHPOfCNT  OF  POSITION  0 
20P  UILL  BE  GREATER  THAN  NUMBER  LISTED',//) 

849  FORNRTdHI,///,  IX, 'DILUTION  OF  PRECISION  -  RCCUNULATIUE  LATITUDE  0 
1ISTRIBUTI0N',/,  IX, 'PROBABILITY  THAT  THE  OOP  UILL  BE  GRERTER  THAN 
2NUMBER  LISTED',//) 


850  FORMATC  1H1, ///,  IX,  'DILUTION  OF  PRECISION  -  ACCUHULATIUE  LATITUDE  D 
1 1STRIBUTIOH1 ,  /,  IX,  'PROBABILITY  THAT  TWEE  D I  MENS  I  OWL  POSITION  OOP 
2  HILL  BE  GREATER  THflN  NUMBER  LISTED',//) 

855  FORHATdHI,///,  IX, 'DILUTION  OF  PRECISION  -  ACCUHULATIUE  LATITUDE  D 
1ISTRI0UTION',/,  IX, 'PROBABILITY  THAT  FOUR  DIIENSIONAL  POSITION  OOP 
2UILL  BE  GREATER  THflH  NUMBER  LISTED',//) 

860  F0RNAT<1H1,///,  IX,  DILUTION  OF  PRECISION  PflARTCTERS  -  ACCUHULATIUE 

1  GLOBAL  DISTRIBUTION*,  //,8X,  NUMBER',  7X,  UOOP',6X,  *  WOP*,  6X, 

2  'WOP',6X,'TDOP',6X,'POOP',6X,'GOOP',/> 

865  FORMATdH  ,9X,F3. 1,3X,6F10.4> 

870  FORHflfT(  1H1,  IX,'  MAXIMUM  fl»  MINIMUM  NUMBERS  SEEN  AT  EACH  LATITUOE 
IS  LONGITUDE’) 

875  FORNAT<  1H0,  /,38 X,  'MINIMUM* ,  //,37X,  'LONGITUDE' ,  //, 

1  32X,2(‘  r,9X),2('2',9X),2(‘3',9X>, /,22X,3(‘5‘,9X,  *0*,9X>, 

2  *5',/,5X,  'LAT',4X,8<'0',QX>, /> 

880  F0RMATdH0,/,38X, 'MAXIMUM', //,37X,  LONGITUDE*,//, 

1  3ZX,2er,9X),2e2',9X),2('3*,9X),/,22X,3('5',9X,'0',9X>, 

2  '5*,/,5X,  ’LAT',4X,8<’0*,flX>,/> 

885  FORMAT< 1H  ,3X,F4.0,3X,36I2> 

800  FORHRT<  1H1,40X,  'NUMBER  OF  SATELLITES' , //,5X,  'LAT',3X,  1617) 

80S  FORNATdH  ,3X,F4.0,5X,  16F7.2) 

900  FORMRTdHO, /,  10X,  'PROBABILITY  (IN  PERCENT)  OF  SEEING  EXACTLY  N  SAT 
1ELLITES',/) 

905  FORMRTdHO,/,  10X, ‘PROBABILITY  (IN  PERCENT)  OF  SEEING  N  OR  MORE  SAT 
1ELLITES',/) 

910  FORMAT(  1H0,/,  IX,  'OH  A  GLOBAL  BASIS  TIE  PROBABILITY  (IN  PERCENT)  TH 
1AT  EXACTLY  N  SATELLITES  HILL  BE  SEEN* ) 

915  FORMAT <  1H0,/,  IX,  ON  A  GLOBAL  BASIS  TIE  PROBABILITY  (IN  PERCENT)  TH 
1AT  N  OR  MORE  SATELLITES  HILL  BE  SEEN*  > 

920  CONTINUE 
END 


SUBROUTINE  Of®  INK  I ) 

C0»10HA)FBIS/P<37,23>/C0h/C<  t8> 

P<I(2WI,2«<1) 

P<I,3HP<I,3X<1> 

P(I(4W<I,4)<C(1) 

P<I,5)"P<I,5X<1) 

P<I,9XIH<P<I,4>> 

P<IJOXOS<P<  1,4)) 

P<l  J1XIN<P<I,3» 

P<l  J2X0S<P<I,3» 

P<  1 ,21  >  -  P<I,22X<6> 

6RA-22808.*C<ft> 

P<l,23)-<SflR/P<l,21)>**1.3 

P<I,6HP<I,21)*<1.-P<I  JHP<I  J» 

P<I,?XQRT<  C<  12>/P<  1 ,6>) 

P<I,8X<12)/<P<I,6HP<I,6)) 

RETURN 

END 


SUBROUTINE  ORBIT<l,T,FER,R,UEL,flC> 

C0NH0N/0RBIS/0P<3?.25)/C0N/C<  18) 

OINENBION  R<3),UEL(3>,RC<3),TRRN8<3,3),Q<3),QEL<3),QC<3),U<36), 

1  E(30),BI0T(36)<LILT<38>JU(36>,CU(38>,SU(36),PER<36> 

RERL  LILT 

IF(T.GT.O)  00  TO  10 
U<IXP<I,3) 

SINE*<SQRT<  1  .-0P<  1 ,  1H*2XIH«JP<  l,5»>/<  1 .  *0P<  1 ,  1  XOS<OP<  1 ,5»> 
RRD-<OP<  1 , 22  >•<  1 . -0P<  I J  W  )  )/<  1 . -KFC I J  XOSCOPC 1 , 3  )  » 
E<IXSIN<SlfC) 

IF<OP<l  J).NE.O.)  00  TO  13 
IF<U<  I  ).LT.<C<3>/2.  ))  E< I  X<  I  ) 

IF«U<I  ).0E.<C<3)/2.  ».RND. 

1  <U<I  >.LE.<3.*C<3>/2.  »)  E<l  X<3HE<I  > 

IFttKI  ).0T.<3.*C<3)/2.  >)  E<IX<4HE<D 
00  TO  20 
13  CONTINUE 

IF<RRO.LT.OP<I,22).RND.U<I).LE.C<3))  E<IX<D 
IF<fiflO.0T.OP<l,22»  E<IK(3K(I> 

IF<RRD.LT.0P<I,22).RND.U<I  >.0T.C<3»  E<IX<4HE<D 
20  CONTINUE 

8I0T<  I  H<PER<  I  >*<E<  I  HCP<  1 , 1 X INE»/<24.*C<4)> 

LILTO  XIGT<I ) 

00  TO  23 
10  0EL>.1 

LILTOXIOTOHT 
E<l  X. 

33  E<l  X<l  HOEL 

V-L I  LT<  I  H<PER<  I  >*<E<  D-OP<  I J  X I  N<E<  I  »  >  >/<24 .  «C<4  >  > 

IF<RB8<V).LE. .00001)  00  TO  30 

IF<V.8T.0)  00  TO  33 

E<l  X<l  HOEL 

OEL-OEL/IO. 

00  TO  35 
30  CONTINUE 
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IF(E<I  >.GE.C(4>>  E<OE<l)-C<4) 

SINU"<SQRT<1.-0P<I  ,1)**2>*SIN<E<I  >>>/<1.-0P<l  <1MC0S<E<I  >>> 
U<l )>ASIN<SINU) 

RP»OP<  1 ,22  M<  1  .-CP<  1 , 1  MCOS<E<  I  >» 

P"0P<I,22M<1.-0P<I  ,1M*2> 

IF<OP<l,  D.NE.O.  )  GO  TO  40 
IF<E(l>.LT.<C<3>/2.»  U<  I  >-U<  I  > 

IF«E<I  ).GE.  <C<3)/2. 

1  <E<I ).LE.<3.*C<3>/2. )))  U(IM(3H(I) 
IF<E<l).GT.<3.<C<3)/2.)>  U<  I  >C<4  MU<  I  > 

GO  TO  25 
40  CONTINUE 

IF<RP.LT.P.RHD.E<I >.LE.C<4)>  U< I  >-U< I > 

IF(HP.GT.P)  IKIK(3H(I) 

IF<RP.LT.P.flW.E<l  >.GT.C<4»  U<  I  >C<4  »U<  I  ) 

23  CONTINUE 

CU<IXOS<U<l)) 

SU<OSIN<U<I  >> 

U<l  MKI  MCP<I,2) 

OP<l  J3M0S<U<I» 

OP<l  ,14>-SIN<U<I» 

SPH-OP<l  ,14W<I,9) 

CPH-SQRT<1.4SPH*SPH) 

OP<l  ,15)-«SIN<SPH>*C<2> 

OP<l  ,18MJ<I) 

F1-1.+0PU  ,1MCU<I) 

Q<1W<I<6VF1 

Q<2M). 

Q<3M). 

OEL<  1  WP<  1 ,1  W<  1 ,7M9U<  I ) 

QEL<2W<I<7MF1 

QEL<3M). 

QC<1)— 0P<I,8MF1*F1 
QC<2)-0. 

QCOMO. 

CflLL  TRNRTX  <TRRNS, I ) 

CflLL  NRfTNUL  <TRRNS,Q,R> 

COLL  NflTHUL  <TRRNS,QEL,UEL> 

CflLL  NRfTNUL  <TRRNS,QC,RC> 

0P<I,24)-<T-  RINT<T+.5»*C<4) 

0P<  1 , 16M1TflN2<R<2  ),  B<  I )  M)P<  1 , 24  ) 

OP<l,  16MP<I,  16MC<2> 

RETURN 

END 


SUBROUT  I*  TRMRTX<TR, I > 

D I HENS I ON  TR<3,3) 

CONHON/ORB I  S/0P<37, 25  >/CON/C<  18) 

TR<1,  IXP<I,  12M0P<I,  13>-0P<I,  11M0P<I,  10MOP<I,  14) 
TR<  1,2  W<  1 , 12  W<  1 , 14  W<  1 , 1 1  W<  1 , 10  W<  1 , 13> 
TR<1<3XP<I<  11M0P<I,9> 

TR<2,1>  0P<I<11W<I  <13W<IJ2W<I  ,10  W<  1 ,14) 
TR<2, 2  >— 0P<  1 ,1 1  W<  I J4  W<  I J2  W<  I  JO  W<  I J3  ) 
TR<2,3>— 0P<  1 , 12  W<  1,9) 

TR<3, 1  W<  1 ,9  W<  1 ,14) 

TR<3,2>  0P<I,9W<I ,13) 

TR<3<3XP<I  JO) 
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TR<1,2>-1.*TR<1,2> 

RETURN 

END 


SUBROUTINE  nRTNUL<T,U,0> 
DIMENSION  T<3,3>,U<3>,0<3> 
00  10  1*1,3 
0(1X1. 

00  10  >1,3 

0(IX<IHT(l,JW<J) 

RETURN 

END 


SUBROUTINE  POINT(RLO,flLfl,TIN,UEC> 
COmON/CON/C(18> 

OirCNSION  UEC<3> 

E1ML0*C(  1  >*C<4)*TIH 
StXLR*C<  1  > 

UEC(  1  X<  10)*C0S<SN)*CQS<EU) 
l€C<2X<  10)*C0S(SN)*SIN<EU> 
UEC(3X(10XIN<SN> 

RETURN 

EI0 


i 

SUBROUT I1C  TflflT  (NRX,MXX,MRfTRIX> 
0IMEN8I0N  NRTRIX< 10020,4) 

DO  10  1*1,10625 

DO  10  11-1,4 

10 

t«ITRIX<l,IIX 

IF(HRX.LT.3>  GO  TO  30 

i 

HXX-<HRX-2  >*(MRX- 1 >*f«X/6 
IttWIRDC-l 

iiivhTU 

MWWftt+1 

NR-0 

DO  25  K*2,MHH 

ICOHC+1 

i 

DO  20  LHCQ,MMN 

KTH.+1 

! 

00  13  1HCT,HN0 

NR-NRH 

HRTRIX<Nfi,2>*K 

MflTRIX(Nfl,3M. 

HRTRIX<Nfi,4M1 

13 

CONTINUE 

1  20 

CONTINUE 

25 

CONTINUE 

SUBROUTINE  UOLUHE  <UUEC,  IO8RT,U0L> 

DIMENSION  UUEC(37,3>,  IOSRT(4),ONE(3>,TUO(3>,T»«EE(3>,FOUR(3> 
OirCNSION  THF(3),THT(3),OHT(3),CROSS(3> 

KR-IDSRTd) 


KB»IDSRT(2) 

KC«IDSRT<3> 

H>)DSAT<4) 

00  10  N-1,3 
0NE(N>UUEC<Kfi,M) 
TM0(NMJUEC(KB,N> 
T>f€E<MMJUEC(KC,M) 
FOUR(NMJUEC(KD,N) 

CALL  UECTOR  (TU0,2,F 
CALL  UECTOR  (THREE, 2 
CALL 

CALL  UECTOR  (TNT,  3,  TIF,  CROSS) 
UOL«ABS<DOT<ONT,  CROSS)) 


SUBROUTINE  COUNRU  <0, IO,NRT,SIG) 
DIMENSION  ID<37),B(4),SIG<6),0<3?,4) 
DIMENSION  RINU(4,4),UKARER(40),TRR(4,4> 
00  20  1-1,4 
00  19  >1,1 
TRR(l,J)-0. 

DO  10  K«1,NHT 
L-ID(K) 

TRR<I,J>TRH(I,J>*G(L,  I  )*0<L,J> 
CONTINUE 

TRA<J,  I  )-TRR<  I ,  J) 

CONTINUE 

TRR<I, I  >-TRR<l,  I  H1.E-12 
CONTINUE 

CALL  LINU2F  (TRA,4,4,RINU, 1,UKARER, IER) 
IFOER.EQ. 120)  00  TO  39 


DO  29  >1,4 
TRR<I,J)-AINU<I,J) 

RINU<l,J>0 

CONTINUE 

CONTINUE 

SIG(1>SQRT(TRR<3,3)) 

SIG<2)-SQRT<TRR<  1, 1HTRR(2,2)) 

S I  OOWWtt  1<SQRT<TAR<  1, 1  »,SQRT<TRfl<2,2)» 
SI6(4>SQRT(TRfi(4,4)) 

SI0<5)-SQRT<TRfl<  1, 1  *TRR<2,2HTAR<3,3» 
8I0(6)-3QRT<TAR<  1, 1HTRR<2,2HTRR<3,3>+TRR<4,4>> 
RETURN 
610(1  >1000. 

SI  0(2  >1000. 

SI  0(3  >1000. 

SI  0(4  >1000. 

SIG(9>1000. 

SI  6(0  >1000. 


comoN/con/c<  is> 

ORTR  C/. 01745329232,57. 295779313,3. 1415926336,6.28318330718 
1 1 . 57079630 , 6076 . 1 16 , 0 ., 0 ., 0 ., 2 . 092ft 1 43504E+O? , 

27. 292 1 1565E-05 , 1.4076380E+ 16, 365. 2563835  ,92.91 
3E+06, 0.0 167272, 23. 44436, -77. 7303, 5290. / 


SUBROUTINE  UECTOROM,  I  ,U2,V3> 
01  MENS  I  ON  U1<3>,U2<3),U3(3> 
00  TO  (10, 15,20),  I 
V3(1MM<1HU2<1> 
U3(2MM(2HU2<2) 
U3<3MM<3HU2<3) 

RETURN 

U3(  1  MI1<  1  >-U2<  1 ) 
U3(2MM(2)-U2<2> 
U3(3)“U1<3>-U2<3> 

RETURN 

U3<  1  MI1(2)*V2<3>-U1<3)*V2<2) 
03(2MM(3>*U2<  1  1  W2<3) 

U3<3>-U1<  1  W2<2MJ1<2W2<  1  > 
RETURN 


FUNCTION  DOT <UI,U2> 

□I HENS I ON  V1<3),U2<3) 

OOT-U 1  <  1  W2<  1  HU  1  <2  W2<2  HU  1  <3  W2<3  > 

RETURN 

END 


SUBROUTINE  UNIUEC  <U,UU) 
DINENSION  U(3),UU<3) 
OENON  -  SORT (DOT  <U,U)) 
UU<1)  -  U(1)/0EN0H 
UU(2)  -  U<2)/DEH0N 
UU(3)  *  U(3)/0EN0N 
RETURN 
END 


THIS  SUBROUTIfC  ALLOWS  THE  USER  TO  SELECTIVELY 
DELETE  UP  TO  TEN  SATELLITES  FROM  THE  NON  INAL  CONSTEL¬ 
LATION.  NOEL1  IS  THE  NUMBER  OF  SATELLITES  TO  BE  DELETED 
AND  IS  ENTEAED  AS  TIC  LAST  I  TEN  ON  THE  FIRST  LINE  OF 
DATA  INPUT.  WEN  SATELLITES  ARE  TO  BE  DELETED,  THE  NEXT 
LINE  OF  DATA  SHOULD  CONTAIN  THE  IDENTIFICATION  NUMBERS 
OF  THE  SATELLITES  TO  BE  DELETED  FROM  THE  CONSTELLATION 
ENTERED  IN  DESCENDING  ORDER  (HIGHEST  ID  NUMBER  TO  LOWEST). 
'JDELd  )*  CONTAINS  TIC  10  NUMBERS  OF  THESE  SATELLITES. 
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ELEVATION  DISTRIBUTION  -  PROBABILITY  THAT  THE  SATELLITE  IN  VIEW  WILL  HAVE  ELEVATION  ANGLES  AS  LISTED 
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DILUTION  OF  PRECISION  -  ACCUMULATIVE  LATITUDE  DISTRIBUTION 
PROBABILITY  THAT  ALTITUDE  DOP  WILL  BE  GREATER  THAN  NUMBER  LISTED 
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DILUTION  OF  PRECISION  -  ACCUMULATIVE  LATITUDE  DISTRIBUTION 

PROBABILITY  THAT  POSITION  DOP  IN  HORIZONTAL  PLANE  WILL  BE  GREATER  THAN  NUMBER  LISTED 
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DILUTION  of  precision  accumulative  latitude  distribution 

PROBABILITY  THAT  THREE  01 MENS I ON AL  POSITION  OOP  WILL  BE  GREATER  THAN  NUMBER  LISTED 


88888! 

090901 


iON«Nin»in-»»NMN«na»-«iiii 
lOaoi'MaistnBnBnfiw-ooi 
I0a«ign'00  0  00g09  09  000i 


19991 

1888! 


oooooooo 


ooooo< 

00000( 

OOOOOI 


!§§^§§§§§8§i 


Isis; 


999900099909990000909990909990099900 


090001 
O  900001 
-  oooooi 


•  g^naor*  —  <oin'»mnr«ca<0  —  oinOQQ9Q< 
iO(o-«*«iMMO»i(*nw-oo909go< 
—  ooooooaooooooc 


o  o 


888888! 

◦  ooooo< 


>04 

is: 


is: 

i  o  c 


ooooo  oo 

8888888 


O 

<N 


in 

a 


O 

n 


in 

n 


o 

* 


ooooooooooo 


8! 
o  < 


>o  < 

18! 


i  O  < 

18! 


OO 

OO 
*■*  r- 
at  (fl 


ooooooooooooooooooooooo 


ooctaMS  n 
--ooo 


m  o  •  r* 

- d 

o 


5588888! 


lO( 

IS! 


1888! 

IOOOI 


>o< 

18! 


ooooooo 

iliiiii 


ooooooooo 

goooooooo 

OOOOOOOOO 

ooooooooo 


oo 

O  00 
O  « 
O  O) 


OOOOOOOOOOOOOOOOOOOOOOOOO 


at  n 
OO 


r*  to 
r*  r* 

d 


Sn  at 
<p  n 


£2g 


oooo 


^  ^  O  < 

o  8  So  < 


§o  4 
O  < 
lOO< 


Hi 


O  O  O  O  < 
0  0  0  0  4 
0  0  0  0  4 


>  O  < 

IS! 


ooooo 

1888888 

lOOOOOO 


OOOOOOO 

OOOOOOO 

ooooooo 

ooooooo 


ooooooo 

000004 
000004 
000004 


oo 


o  o 

8? 


OOOOOOOOOOOOOOOOOOOOOOOOO 


O  01 
in  d 
a>  r- 


at  ao 

o 


o  oo  n 
in  ♦ 
o  o  o 


in<0oco«*-aa» 

oooooooo 


888888 


Igoo 


>09 


*8! 


!  111! 


>0  0  0004 
1000004 
>000004 


UJ  O 

a  m 

o 


<  in 
-j  in 


oooooooooo ooooo ooooooooooooooooooooo 

22S2229229hfllNW"#,flop‘w®*p*,fl#N«,0,",nn’"OOOQ 
02222222920,toO«doo«d<bin<r*4*jdd*-**  —  6000000 
OOOOOOOOO  OO»0Q(Onoi«-«- ooooo  oooooooooooooo 

-  ~  ~  -  . oooo 


00000000004 

00000QOO004 

00000000004 


iiiiiiiiiiiiiiiiii 


o 

<0 


in 

(O 


o 

r* 


tn 

r* 


o 

ao 


w 

as 


o 

at 


oooooooooooooooooooooooooooooooooooo 

2922222222?,rr®5"®"«w00'D*w9|flnNN’,,QOOpo 
ogogooOOOO4ni0Ot^O«4»>-r»**-t-.*d~SQQQOOOQQQQ 
000  0000000««KI(inwwN*-*-*'00  0  00000000000  0 

OOOOO 


>0004 

10004 

10004 


060006060066666006006000606600666600 

f2222  22Q222^nOM0lOl0«^NOn«(1dln^lQOOooooo 
00600009  OQ^r4D4Ddaj^aonaiin»-r-^c<4gQo5g6Q65  0Q 
OOOOOOOOOO4BI0iAtnnNM^*«*OOOOOOOOOOOOOO 


ooooooo 


!!!§! 


>04 

18! 


00000000000000000000000000000000000 

§S2222222299©*~wp*®r*,c0*~,®>r,",,2922O2O999 
222222299o9P*f',fswf'4*"0®r'*incjM~gggggggggo 
600000d0000«t^'----0000000000000000 


>04 

*  O  4 
)  O  4 


ooooooooo 

Slliiililii 


999999990099900990999999090009999999 

OOQOQOOPOOOOainaieBear'(ainnr«''OOQQ0099000QO 

2S92952999QOK>M^o9«o»'nt«oo5ooooogooo 

OOooo99999ooaifltff»nn«»'Oooo999oooo9ooo 

------------  999999999999 


X  ONt<oa9N«Ba9N*oa9nttfaON«ioa9Nvioa9Nfils9 
Z  9  --. . 


'NNH««nnnnni 


iinmuitnuxoisvior- 


0  000  0.000  0 


DILUTION  OF  PRECISION  -  ACCUMULATIVE  LATITUDE  DISTRIBUTION 
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PROBABILITY  (IN  PERCENT  I  OF  SEEING  EXACTLY  N  SATELLITES 
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tput  -  Satellite  User 
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Appendix  E 

Selected  Data  Extracts 


This  section  of  the  appendix  contains  some  selected  output  from  the  numerous  computer  runs 
made  and  used  in  the  analysis  of  the  baseline  constellation  and  its  modification.  The  "*nplete 
program  output  listing  is  provided  for  the  global  distribution  run  made  for  each  of  these  satellite 
constellations.  Portions  of  the  output  from  several  otter  runs,  which  were  used  in  analyzing  the 
effect  of  satellite  losses  on  each  constellation,  have  also  been  provided.  Finally,  extracts  of  output 
from  several  of  the  computer  runs  used  in  the  analysis  of  the  geometric  performance  for  the 
space-based  user  are  included  as  indicative  of  the  performance  obtainable  for  varying  antenna 
half-angle  beamwidths. 
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DILUTION  OF  PRECISION  -  ACCUMULATIVE  LATITUDE  DISTRIBUTION 
PROBABILITY  THAT  ALTITUDE  OOP  WILL  BE  GREATER  THAN  NUMBER  LISTED 
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THE  SATELLITE  MOST  NEARLY  OVERHEAD  IS  USED  AS  ONE  OF  THE  FOUR 
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